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PREFACF, 


This  is  the  Standard  Ocean  Evaluation  Group's  final  report  on 
Its  search  for  a  data  retrieval  system  that  will  be  installed  as 
Standard  Ocean  in  the  data  bank  of  the  Long  Range  Acoustic  Propaga¬ 
tion  Project  (LRAPP),  sponsored  by  the  Naval  Ocean  Research  and  De¬ 
velopment  Activity  (NORDA).  LRAPP  formed  the  evaluation  group  in 
January  1979  to  assemble  a  roster  of  candidate  systems  and  evaluate 
each  against  project  requirements.  The  group  comprised  J.  G.  Colborn, 
chairman  (Naval  Ocean  Systems  Center);  S.  C.  Daubin,  Jr.  (Pacific- 
Sierra  Research  Corporation  ( PSR  ] ) ;  E.  liashimoto  (NORDA);  and  F.  J. 

Ryan  (Ocean  Data  Systems,  Inc.  [ ODS I ] ) .  Standard  Ocean  is  to  be 
installed  early  in  FY  1981. 

Preliminary  results  of  the  evaluation  group's  search  were  reported 
in  a  memorandum  to  I.cdr.  Kirk  Evans  (NORDA)  and  John  H.  Locklin  (ODSI) 

it 

on  JO  April  19/9.  The  present  report  reflects  subsequent  discussions 
held  at  NORDA,  further  evolution  of  the  prime  candidate  model,  CDEM, 
and  recent  evaluation  efforts. 

The  draft  of  the  report  was  submitted  for  review  in  August  1979 
(results  are  current  as  of  that  date);  the  report  was  approved  for 
distribution  on  7  February  1980.  In  the  interim,  the  following  changes 
occurred:  LRAPP  was  renamed  the  Surveillance  Environmental  Acoustic 

Support  (SEAS)  Project;  and  Fleet  Numerical  Weather  Central  (FNWC)  is 
now  called  the  Fleet  Numerical  Oceanography  Center  (FNOC).  Readers 
will  note  that  the  text  retains  the  former  nomenc lature .  Also  in  the 
meantime,  F.  !.  Kvatt  changed  his  affiliation  to  Science  Applications, 
t  nc  . 

The  efforts  of  the  following  individuals  at  PSR  contributed  to 
v  h  i  s  report:  Christine  D’Arc  edited  it;  Joan  Pederson  typed  it;  Laurie 
biackeby  arid  Timothy  Had  lock  prepared  the  artwork. 

* 

"Prc 1 iminary  Evaluation  of  Candidates  for  the  Standard  Ocean 
Retrieval  System  and  Next  Steps  Toward  implementation." 
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SUMMARY 


This  report  describes  and  evaluates  eight  existing  or  proposed 
oceanographic  models  as  candidates  for  Standard  Ocean,  a  data  retrieval 
system  to  be  installed  in  the  Long  Range  Acoustic  Propagation  Project 
(LRAPP)  data  bank.  The  primary  purpose  of  Standard  Ocean  is  to  pro¬ 
vide  range-dependent  sound-speed  profiles  for  input  to  NORDA's  numerical 
acoustic  models.  Standard  Ocean  will  also  be  used  to  support  the  ob¬ 
jective  analysis  of  environmental  data  collected  during  exercises  at 
sea.  The  candidate  systems  and  their  parent  organizations  are  as 
follows : 


AUTO-OCEAN  (NORDA) 
FIB/EOTS/EXTRA  (FNWC) 

C.DEM  (NAVOCEANO) 

ICAPS  (NAVOCEANO) 

ODSI  (Ocean  Data  Systems,  Inc.) 


CFDI.  (Geophysical  Fluid  Dy¬ 
namics  Laboratory,  NOAA, 
Princeton  University) 

HYDAT  (FNWC) 

SIMAS  (NUSC/New  London) 


The  Standard  Ocean  Evaluation  Group  assessed  each  candidate  accord¬ 
ing  to  criteria  indicated  in  the  following  description  of  desired 
Standard  Ocean  capabilities.  Standard  Ocean  is  to  provide  accurate, 
realistic,  and  seasonal  (preferably  monthly)  surface-to-bottoa  profiles 
of  sound  speed,  temperature,  and  salinity  it*  each  oceanic  1*  x  1* 
square.  The  sound-speed  profiles  should  be  in  a  format  suitable  for 
numerical  acoustic  models.  The  profiles  sh Mid  accurately  reproduce 
all  acoustically  significant  features.  The  degree  of  oceanic  vari¬ 
ability  in  each  square  should  be  indicated.  Standard  Ocean  should 
operate  rapidly  and  inexpensive  v.  it  should  be  easily  usable  by  the 
nonspecialist.  Finall*',  the  candidate  chosen  should  ie  competitive 
in  acquisition  cost  and  availability. 

Using  those  criteria,  we  found  the  following  six  candidates  un¬ 
suitable  for  Standard  Ocean:  FIB/EOTS/EXTRA,  GFDL,  HYDAT,  ICAPS,  ODSI, 
ami  SIMAS.  The  remaining  two,  AUTO-OCEAN  and  CDKM,  met  or  exceeded  most 
Standard  Ocean  criteria.  Nevertheless,  each  requires  modification 
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bcfore  it  can  be  adoptee  as  Standard  Ocean-  AUTO-OCEAN,  an  operational 
retrieval  system,  requires  a  finer  grid  spacing  for  its  sound-speed 
profiles  and  the  addition  of  temperature  and  salinity  profiles  as  out¬ 
puts.  GDEM,  an  objective  analysis  model  being  developed,  requires 
refinement  cf  its  objective  analysis  technique  to  remove  anomalies  in 
the  middepch  sound-speed  profiles.  The  developer  is  redesigning  part 
oi  GDEM  to  correct  those  faults.  Completion  of  the  GDEM  analysis  for 
the  North  Atlantic,  North  Pacific,  and  Mediterranean  is  expected  by 
early  FY  1981.  The  Indian  Ocean  portion  is  to  be  completed  at  an  un¬ 
specified  future  time.  Should  GDEM  not  be  ready  by  early  FY  1981,  the 
target  date  for  installation  of  Standard  Ocean,  we  recommend  that  AUTO¬ 
OCEAN  be  used  in  the  interim,  revised  in  key  areas  of  LRAPP  interest 
as  indicated  above.  Even  when  GDEM  is  completed,  we  recommend  storing 
the  GDEM  outputs  in  the  AUTO-OCEAN  framework  to  satisfy  the  require¬ 
ment  that  acoustic  model  inputs  be  automated. 

Me  candidate  has  the  required  capability  of  indicating  oceanic 
variability.  We  recommend  that  LRAPP  investigate  how  that  capability 
could  be  developed  for  Standard  Ocean  after  installation. 
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I .  INTRODUCTION 


This  report  presents  the  results  of  a  search  for  a  data  retrieval 
system  to  be  installed  as  the  Standard  Ocean  in  LRAPP's  data  bank. 

The  Standard  Ocean  Evaluation  Group  determined  that  eight  systems 
warranted  consideration  as  Standard  Ocean  candidates,  and  wt  examined 
each  in  light  of  LRAPP's  requirements. 

THE  NEED  FOR  STANDARD  OCEAN 

Standard  Ocean  was  conceived  by  Lcdr.  Kirk  Evans  (NORDA)  to  meet 
LRAPP's  need  for  the  rapid  retrieval  of  accurate,  realistic  range- 
dependent  oceanographic  data,  primarily  sound-speed  data,  for  input 
in  NORDA 's  numerical  acoustic  models.  LRAPP  uses  those  models  for  a 
variety  of  purposes  including  exercise  preassessment,  exercise  post¬ 
analysis,  and  area  assessments.  Standard  Ocean  is  also  needed  to 

provide  inputs  for  the  objective  analysis  of  exercise  environmental 

* 

data.  Other  potential  users  of  Standard  Ocean  include  Fleet  Numerical 
Weather  Central  (FNVC),  NORDA  Code  320,  and  the  oceanographic  community 
at  large. 

Standard  Ocean  will  furnish  a  variety  of  oceanographic  data  besides 
sound  speed.  It  should  be  able  to  provide  accurate,  realistic  surface- 
to-bottom  profiles  of  sound  speed,  temperature,  and  salinity  in  each 
1*  *  1‘  square  of  the  oceanic  northern  hemisphere.  Each  profile  will 
be  typical  for  a  particular  square  in  one  of  the  four  seasons  (if  attain¬ 
able,  monthly  resolution  is  desired).  Tike  profiles  will  represent  the 
most  probable  oceanographic  conditions  in  the  given  square  and  time 
period.  "Most  prcoable"  does  not  signify  average  or  mean;  though  sta¬ 
tistically  correct,  an  average  profile  might  never  be  observed  in  that 
location. 

Standard  Ocean  will  produce  acoustically  <*  .urate  sound-speed 
profiles.  Hi  at  is,  sound-speed  magnitudes,  channel  axes  and  layer 

V ' 

J .  Locklin  et  al.,  Ui AFP  Jbj  eo  i  :  oe  A  AVrteftJ  of  Ob¬ 
jective  Ar+ilysm  :k'htOKed  for  the  LRAPP  LktUt  Prix/risn,  LRAFF. 

20  March  197V. 
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depths,  and  horizontal  and  vertical  gradients  should  be  those  most 
likely  to  be  observed  at  the  location  and  time  specified.  The  sound- 
speed  profiles  should  be  in  a  format  suitable  for  numerical  acoustic 
models . 

Though  not  used  directly  in  numerical  acoustic  models,  temperc * 

Cure  and  salinity  data  are  valuable  for  verifying  the  quality  of  c >  ' 
output  that  ia  used.  The  accuracy  of  sound-speed  output  cannot  be 
determined  from  sound  speed  alone.  Other  more  stable  measures  of 
accuracy  are  needed.  Temperature  and  salinity  change  in  highly 
dicrable  ways  (i.e.,  away  from  the  oceanic  boundaries,  heat  an;  salt 
are  not  added  or  subtracted  from  the  ocean  except  in  the  s  'face  layer 
above  100  m,  in  some  shallow  areas,  and  to  a  much  lesser  extent  on  the 
ocean  bottom).  Therefore  temperature  and  salinity  data  tacilitate 
quality  checks  of  model  output,  two  common  ones  being  (1)  that  density 

does  not  decrease  with  depth,  and  (2)  that  temperature  versus  salinity 

* 

versus  depth  follow  known  empirical  relationships.  Hence  the  need 
for  temperature  and  salinity  data  in  Standard  Ocean. 

The  degree  jf  oceanic  variability  in  each  square,  though  not  of 
immediate  concern  to  acoustic  modelers,  should  eventually  be  added  to 
Standard  Ocean's  capabilities.  Over  long  periods,  variability  in  the 
upper  kilometer  of  the  ocean  can  eaust  marked  differences  between  pre¬ 
dicted  and  observed  conditions.  For  Instance,  the  variability-caused 
"noise"  in  tiie  observed  surface  temperature  can  be  as  great  in  ampli¬ 
tude  as  tiie  annual  signal.  Variability  outputs  needed  for  Standard 
Ocean  include  measures  of  the  incidence  of  sonic  layers,  their  minimum 
and  maximum  depths,  the  incidence  of  eddies  and  fronts,  and  the  minimum 
and  maximum  ranges  of  temperature,  salinity,  and  sound-speed  profile*-.. 

Standard  Ocean  vill  be  an  automated  system  that  a  nonoccanograniier/ 
notueousi ielan  can  use  to  quickly  and  cheaply  select  the  appropriate 
sound-speed  prufiies  for  model  runs.  in  effect.  Standard  Ocean  will 

The  rationale  lot  temperature  and  salinity  analysis,  and  the  pro¬ 
cedures  involved,  are  explained  ‘a  greater  detail  in  o.  I.  Mamayev, 
."ejrTj'cnKurc-.'i.!.'  mi  t:(  Atulua ia  •.>/*  i/cr.’J  »u f *r rv ,  Elsevier  Ocean¬ 

ography  Serlea,  Vol .  11,  Elsevier  Press,  Amsterdam,  The  Netherlands, 

19  7^,  and  in  Appendix  F. 
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remove  the  oceanographer  from  the  data  selection  "loop"  and  replace 
him  with  a  computerized  data  retrieval  system.  The  oceanographer’s 
attention  can  thus  be  freed  for  model  'runs  requiring  more  specialized 
inputs  and  interpretations. 

Standard  Ocean  will  function  as  a  separate  unit  within  the  LRAPP 
data  bank.  Figure  1  depicts  Standard  Ocean's  role  in  the  flow  of 
information  from  data  sources  to  acoustic  models.  Environmental  data 
from  sources  such  as  NODC  and  FNWC  are  screened  and  edited.  Depending 
on  the  technique  employed,  either  the  most  probable  real  profiles  are 
select ad,  or  all  observed  data  are  objectively  analyzed  to  produce  a 
smoothed,  representative  field.  The  analyzed  data  become  the  Standard 
Ocean.  Next,  the  data  are  retrieved  from  Standard  Ocean  and  combined 
with  other  environmental  data  such  as  shipping  distributions  and  bathy¬ 
metry  for  input  first  into  the  appropriate  reformatting  software  and 
then  into  acoustic  models  such  as  PE,  FACT,  and  ASTRAL.  Alternatively, 
Standard  Ocean  can  input  data  into  the  objective  analysis  module  of  the 
data  bank  and  then  into  the  reformatting  modules  and  acoustic  models. 

In  that  capacity,  Standard  Ocean  provides  the  first-guess  climatology 
necessary  to  initiate  the  objective  analysis  of  a  particular  data  set. 

STANDARD  OCEAN  REQUIREMENTS 

The  foregoing  description  of  the  attributes  desired  in  Standard 
Ocean  yielded  the  following  list  of  requirements,  which  we  used  to 
evaluate  the  candidate  systems: 

•  Output  products  and  technical  quality .  Standard  Ocean  should 
provide  accurate,  realistic,  and  seasonal  (preferably  monthly) 
surface-to-bottom  profiles  of  the  most  probable  sound  speed, 
temperature,  and  salinity  in  each  oceanic  1°  x  1°  square. 

The  sound-speed  profiles  should  be  in  a  format  suitable  for 
numerical  acoustic  models  and  should  accurately  reproduce  all 
acoustically  significant  features.  The  degree  of  oceanic 
variability  in  each  square  should  be  indicated  whenever  pos¬ 
sible.  Standard  Ocean  should  also  produce  outputs  that  can 
be  used  to  support  the  objective  analysis  of  synoptic  exercise 
environmental  data. 


Standard  Ocean  in  LRAPP  data  bank 


•  Level  of  automation.  Standard  Ocean's  output  should  be  easily 
usable  by  a  nonspecialist. 

•  Turnaround  time  and  operating  cost.  Standard  Ocean  should 
have  minimal  turnaround  time  and  operating  cost. 

•  Time  and  ooet  to  acquire.  Standard  Ocean  should  be  competitive 
in  time  and  cost  required  to  incorporate  in  the  LRAPP  data  bank. 

Section  II  presents  the  evaluation,  which  concludes  with  our  find¬ 
ings  on  the  candidates  judged  unsuitable  or  suitable  for  Standard 
Ocean.  In  Sec.  Ill,  we  narrow  the  choice  further  and  consider  the  tasks 
that  remain  before  Standard  Ocean  can  be  installed. 
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II.  EVALUATION 

Tills  section  lists  the  eight  Standard  Ocean  candidates  we  evalu¬ 
ated,  describes  their  distinguishing  characteristics,  then  explains 
why  we  find  certain  candidates  unsuitable  and  others  suitable  for 
Standard  Ocean. 

THE  CANDIDATES 

The  following  list  names  the  parent  agency,  key  persons  involved 
in  each  model's  development  and  operation,  and  the  appendix  in  this 
report  that  describes  the  system  more  fully. 

•  AUTO-OCEAN,  in  operation  at  NORDA,  Code  320,  NSTL  Station,  Bay 
St.  Louis,  Mississippi  (Appendix  A). 

•  FIB/ EOTS / EXTRA  (Fields  by  information  Blending/Extended  Ocean 
Thermal  Structure/EXTRAction) ,  developed  by  Manfred  Holl  of 
Meteorology  International,  Inc,,  for  FNWC,  Monterey,  Cali- 
for-  a  (Appendix  B) . 

•  GDF.M  (Generalized  Digital  Environmental  Model),  developed  by 
Thomas  M.  Davis  at  NAV0CEAN0,  NSTL  Station,  Bay  St.  Louis, 
Mississippi  (Appendix  C) . 

•  GFDL  (ue  model  name  given),  developed  by  Sydney  Lev  It  us  and 
Abralau.  It.  Oort  of  the  Geophysical  Fluid  Dynamics  Laboratory, 
NO.AA  at  Frinceto.  University,  Princeton,  New  Jersey  (Appendix 
D). 

•  HYDAT  (tiYdro-Ci  imatological  DATs  Basel,  developed  by  Ocean 
Data  Systems,  It.o.,  Montetey,  California,  for  FNWC;  operated 
and  taodiJied  b,  ’.velyn  A.  Hes_>,  FNWC  (Appendix  E)  . 

•  U'AFS  (integrated  Coasaand  Antisubmarine  Warfare  I  red i c t ion 
System',  developed  by  Alvan  Fisher,  Jr.  and  operated  by  William 
R.  F'oyd  and  Haul  Moersdorf  of  NAVOCEANO  (Appendix  F) . 

•  ODSi  (no  model  name  g'vcu),  proposed  model  to  be  developed 
by  Capt .  Paul  M.  Wolft  o!  Ocean  Data  Systems,  Inc.  (Appen¬ 
dix  G). 
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•  SIMAS  (Sonar  In-Situ  Mode  Assessment;  System),  developed  by 
Eugene  M.  Podeszwa  of  NUSC,  New  London,  Connecticut  (Appen¬ 
dix  H)  . 


Distinguishing  Characteristics 

The  candidates  are  differentiated  by  several  key  characteristics, 

summarized  in  Table  1.  One  is  the  method  used  to  derive  the  model's 

output  products.  F IB/ EOTS/ EXTRA ,  GDEM,  GFDL,  and  ICAPS  use  an  objec- 

* 

tive  analysis  technique.  AUTO-OCEAN,  HYDAT,  and  SIMAS  store  and  re¬ 
trieve  typical  observed  prof  * \es  for  each  location  and  time  period;  in 
addition,  AUTO-OCEAN  performs  some  smoothing  of  typical  profiles. 

ODSI's  analysis  technique  is  unclear. 

The  objective  analysis  approach  has  the  advantage  of  working  well 
in  areas  of  sparse  data,  provided  the  surrounding  areas  have  adequate 
data.  Difficulty  may  arise  in  analyzing  areas  where  all  the  data  were 
collected  over  a  short  period  or  only  during  one  or  two  seasons 
(seasonal  aliasing).  The  objective  analysis  technique  can  also  falter 
in  analyzing  situations  where  two  or  more  water  masses  occupy  a  single 
square.  If  the  masses  have  very  distinct  characteristics,  the  analysis 
may  fall  Into  the  trap  of  producing  an  average  result,  which  may  never 
be  found  in  reality. 

The  typical  observed  profile  approach  works  well  in  areas  with 

abundant  observations,  yhere  it  produces  representative  profiles  close 

to  the  most  probable  real  conditions.  However,  fewer  than  3  million 

t 

deep-ocean  measurements  have  been  taken  worldwide,  and  they  cluster 
(1)  along  shipping  lanes  and  around  ocean  weather  stations  in  the 
northern  hemisphere ,  and  (2)  in  the  warmer,  less  stormy  months.  When 
profiles  are  desired  for  the  times  and  places  unrepresented  by  data, 

* 

Ttie  objective  analysis  models  vary  in  the  sophistication  and  range 
of  their  techniques.  The  common  thread  is  that  their  outputs  are  not 
observed  profiles  but  profiles;  statistically  derived  from  real  data. 
Typically  they  are  smoothed  approximations  of  average  or  most  probable 
oceanographic  conditions  and  do  not  contain  the  "noise"  found  in  real 
data.  Lock  1  in  et  al.  review  the  objective  analysis  techniques  of  GDEM, 
FIB/EOTS/EXTRA,  and  GKDL. 

•+* 

As  opposed  to  surface  measurements,  of  which  there  are  about 
10  allilon. 
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Table  1 

DISTINGUISHING  CHARACTERISTICS  OF  CANDIDATE  MODELS 


Cl  laracter  istlc  Model 

Method  of  deriving  outputs: 

Objective  analysis  .  FI B/EOTS /EXTRA,  GDEM, 

GFDL,  I CAPS3 

Typical  observed  profiles  .  AUTO-OCEAN,  HYDAT, 

SIMAS 

Number  of  layers  represented: 

One  .  HYDAT,  ODSI 

More  than  one  .  AUTO-OCEAN,  FIB/EOTS/ 

EXTRA,  GDEM,  GFDL , c 
ICAPS,  SIMAS 

FIB/EOTS/EXTRA,  GDEM 
(temperature  and 
salinity  only) ,  GFDL 
AUTO-OCEAN,  GDEM  (sound 
speed  only),  HYDAT, 
ICAPS,  SIMAS,  ODSI 

q 

ICAPS  uses  a  hybrid  approach,  combining  XBT  observations 
with  an  analyzed  temperature  profile  after  identifying  the 
local  water  mass. 

^Hand-smoothed . 
c 

GFDL  outputs  32  horizontally  analyzed  surfaces  of  temper¬ 
ature,  salinity,  sigma-T,  and  oxygen. 


Greatest  depth  of  output  products: 
<  5000  m  . 

Bottom  . 


the  skills  of  the  oceanographer  are  strained  and  he  may  be  forced  to 
use  "artistic  license"  in  place  of  scientific  techniques.  Assuming 
that  real  data  offer  the  most  accurate  standard,  we  evaluated  the 
model  outputs  of  GDEM  by  comparing  them  with  typical  observed  profiles, 
particularly  at  locations  where  the  data  are  plentiful  (see  Appendix  C) . 
Yet,  to  be  effective  worldwide,  Standard  Ocean  may  need  to  incorporate 
an  objective  analysis  technique  for  use  in  sparsely  sampled  areas. 

The  second  distinguishing  characteristic  is  the  number  of  layers 
used  to  represent  variables  such  as  temperature  and  salinity.  ODSI 
and  HYDAT  use  a  single-layer  model.  The  rest  of  the  candidates  are 
multilayer.  The  multilayer  approach  assumes  that  the  water  in  the 
upper  layers  (typically  <400  m)  differs  from  deeper  water  in  dynamic 
characteristics  and  time  constants.  AUTO-OCEAN,  GDEM,  FIB/EOTS/EXTRA, 
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SIMAS,  and  I CAPS  employ  a  two-  or  three-layer  model.  GFDL  outputs 
horizontally  analyzed  surfaces  at  32  depths  of  temperature,  salinity, 
sigma-T,  and  oxygen.  Unlike  the  other  multilayer  models,  GFDL  does  not 
perform  a  vertical  analysis,  although  the  values  at  the  different  depths 
can  be  recombined  to  form  vertical  profiles.  Using  a  multilayer  model, 
the  modeler  can  simultaneously  represent  the  deeper  ocean  with  an 
annual  or  w inter-suraner  resolution  and  the  upper  layers  with  seasonal 
or  monthly  resolution.  (ICAPS  is  an  exception  in  this  respect.  It 
merges  expendable  bathythermograph  (XBT)  data  with  an  appropriate  deep 
(>200  ra)  temperature  profile  after  determining  the  local  water-mass 
characteristics — see  Appendix  F. )  In  the  final  output  products  the 
layers  are  merged  over  a  range  of  depths  around  their  overlap  or  meet¬ 
ing  point.  Another  reason  for  the  importance  of  temperature  and  salinity 
in  Standard  Ocean  (besides  the  reasons  given  in  Sec.  I)  is  that  errors 
are  easier  to  detect  in  merges  based  on  temperature  and  salinity  than 
in  merges  based  on  sound  speed  alone. 

The  third  distinguishing  characteristic  is  the  greatest  depth  to 
which  the  analysis  extends.  FIB/EOTS/EXTRA  and  GFDL  produce  output 
products  no  deeper  than  5000  ra.  The  rest  extend  to  the  bottom,  except 
GDEM’s  temperature  and  salinity  models,  which  do  not  extend  beyond 
800  in.  The  depth  limitation  does  not  exclude  FIB/EOTS/EXTRA  and  GFDL 
from  consideration,  since  their  output  could  be  extrapolated  to  the 

it 

bottom  with  existing  deep  sound-speed  data. 

Output  Products 

Table  2  summarizes  the  output  products  and  coverage  of  the  candi¬ 
date  models.  The  detail  of  the  descriptions  corresponds  with  the 
amount  of  information  available.  For  example,  much  more  information 
is  available  about  GDEM  chan  about  HYDAT. 


evaluation 

Our  evaluation  of  tiie  candidate  models  against  the  requirements 


J.  J .  Audet  ,  Jr..  XAtain.il  cVt' t  u\ 
Channel  In  t  cr/Vrcmv  for  the  'Wan  and  . 
ay  he  re ,  Acoustic  Environmental  Support 
ideal  Note  TN-Z'i-Oi,  April  1975. 


li  ‘  V;  th  i'ftartr  an. I  i>eey  Sound 
Sea  a  of  the  Northern  Henri - 
Detachment,  ONR,  AESD  Tech- 


OUTPUT  PRODUCTS  AND  COVERAGE  OF  CANDIDATE  MODELS 
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speclfled  In  Sec.  I  revealed  that  each  candidate  has  unique  capabilities 
meeting  one  or  more  of  LRAPP's  requirements;  none  meets  all  of  them. 

Unsuitable  Candidates 

F IB /EOTS /EXTRA  is  a  documented  system  both  undergoing  testing  and 

in  operation  at  FNWC.  It  is  used  on  a  CDC  6500  to  objectively  analyze 

oceanic  thermal  structure  with  XBT  observations  as  inputs.  The  outputs 

are  "fields, "  or  horizontally  analyzed  surfaces  of  temperature  and  its 

so-called  first  and  second  differences  defined  on  FNWC's  standard 

2  2 

63  x  63  grid  covering  the  northern  hemisphere  (i.e.,  AT/AZ  and  A  T/AZ  ). 
The  fields  are  computed  at  26  depths  between  the  surface  and  5000  m, 
with  increased  resolution  near  the  bottom  of  the  mixed  layer.  Spatial 
resolution  of  the  63  *  63  grid  at  60*N  is  ~  380  km.  For  limited  geo¬ 
graphic  areas,  a  grid  of  higher  resolution  is  available  down  to  1/8  the 
standard.  The  system  accepts  as  inputs  temperature  profile  .  that  are 
irregularly  spaced  in  time  and  space.  Sound  speed  is  computed  usin,u 
analyzed  temperature,  archived  salinity,  and  the  FNWC  program  EXTRA. 

In  regions  lacking  XBT  observations,  the  model  reverts  to  the  archived 
climatology . 

The  disadvantages  of  FIB/ EOTS /EXTRA  are  that  (I)  it  analyzes  tem¬ 
perature  only,  using  Indirect  means  to  compute  sound  speed,  (2)  the 

analysis  is  sensitive  to  the  sparsity  of  synoptic  XBT  observations,  so 

* 

it  relies  heavily  on  static  climatology,  (3)  the  merge  of  the  surface 
temperature  analysis  (0-400  m)  with  the  deep  climatology  could  produce 
unrea 1 ist le  results,  and  (4)  the  analysis  extends  onlv  to  5000  m  (not 
a  serious  deficiency). 

The  advantage  of  FIB/ EOTS /EXTRA  is  the  ease  with  which  its  objec¬ 
tive-  analysis  technique  could  be  developed  for  Standard  Ocean  by  (1) 
improving  the  spatial  grid  to  at  least  1*  *  1*,  (?)  increasing  the 
depth  coverage  to  possibly  11,000  ®,  and  (3)  improving  the  input  cli¬ 
matology  and  adding  salinity  as  an  output.  It  would  be  most  effective 
to  run  FIB /EOTS/ EXTRA  at  FNWC  and  transfer  the  outputs  to  the  !iA?f 


* 

This  problem  could  be  alleviated  by  (1)  increasing  the  observa¬ 
tion  period  fioo  the  present  12  hr  to  approximately  one  month  and  (2) 


updating  the  climatology  with 
Cion  Data  Set  (MOODS)  contain 


FNWC's  new  Master  Oceanographic 
icg  ~  2.8  x  106  observations. 


Observa- 
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d.ira  bank  for  inclusion  in  Standard  Ocean.  The  likely  cost  of  such 
implementation  is  unknown. 

('.F'nh  is  an  objective-analysis  model  that  produces  horizontal  fields 
of  analyzed  temperature,  salinity,  sigma-T,  and  oxygen  in  worldwide 
la  '  1°  squares.  Temporal  resolution  Is  monthly.  The  data  are  output 
at  12  depths  from  the  surface  to  5000  m.  GFDi-’s  most  serf our  flaw  is 
its  lack  of  attempt  to  preserve  the  vertical  integrity  of  the  output 
fields.  Sound  speed,  not  a  direct  output  product,  must  be  drived  from 
the  model's  temperature  and  salinity  outpucs.  As  vertical  integrity  is 
not  preserved  in  those  parameters,  the  sonic  layer  deptns  and  secondary 
channels  in  any  derived  sound-speed  profile  will  be  distorted.  In  ad¬ 
dition,  the  model  produces  seasonal  aliasing  in  squares  lacking  cora- 
ple’r  seasonal  data.  It  would  probably  be  difficult  to  modify  the  model 
for  Standard  Ocean  application  because  the  objective-analysis  computer 
code  is  specific  to  the  GFDL  computer;  the  output  data  format  may  also 
be  unique. 

The  main  advantages  of  GFDL  are  that  the  initial  analysis  of  the 
fields  is  complete,  the  coverage  is  worldwide,  and  the  model  is  pre¬ 
sumably  available  now. 

•<’  .  i  s  a  retrieval  system  for  selecting  the  most  representat !' 

observed  temperature  and  salinity  profiles  In  a  given  area  (multiple 
of  l  ‘l  squares).  Coverage  is  worldwide  at  intervals  of  one  or  two 
months.  [’rutiles  extend  f  rom  surface  to  bottom. 

HYDAT  is  unsuitable  because  (l)  its  analysis  does  not  attempt  to 
preserve  representative  sonic  layer  depths,  (2)  it  sometimes  fails  to 
d i f t erett t  la t e  between  two  water  masses  (e.g.,  when  water  masses  arc 
different  at  the  surface  but  the  same  at  241  m,  and  when  water-mass 
characteristics  reverse  between  the  surface  and  243  a),  and  (3)  it 
does  not  preserve  the  horizontal  integrity  of  the  output  fields. 

HYDAT's  older  data  set  of  750,000  observations  Is  to  be  replaced 
in  late  19 7v  by  FNWC'a  new  MOODS  observation  data  set.  Evelyn  A.  Hess, 
who  is  debugging  and  modifying  HYDAT  at  FNWC.  does  not  consider  the 
model  ready  for  release  yet.  As  HY'DAT  is  designed  to  operate  with 
F5iUC  software,  it  may  requite  modification  to  operate  on  the  LKAFF  data 
bank.  finally,  an  experienced  operator  is  needed  to  run  HY'DAT. 

The  main  advantages  of  HY'DAt  are  its  worldwide  coverage  and  its 
ability  to  sometimes  identify  tlie  presence  of  two  water  masses  in  an 
area  and  give  their  mean  character  1st  lea . 
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ICAPS  Is  a  shipboard  system  for  predicting  sonar  performance  based 
on  XBT  inputs,  archived  water-mass  characteristics,  and  FACT  model  runs. 
ICAPS  merges  an  observed  XBT  with  an  appropriate  deep  (>200  m)  tempera¬ 
ture  profile  by  identifying  the  water  mass  in  the  upper  200  m.  The 
temperature  and  annual  salinity  profiles  are  used  to  compute  a  sound- 
speed  profile,  which  is  fed  into  the  FACT  model  to  produce  acoustic 
eigen-rays  and  propagation  loss. 

ICAPS  is  operational  only  in  the  northern  hemisphere.  It  cannot 
predict  what  the  range-dependent  environmental  field  will  be  during  a 
particular  period;  it  requires  real  XBT  data  as  inputs.  Furthermore, 
it  assumes  a  uniform  ocean  in  the  vicinity  of  the  XBT  observation. 

As  a  result,  its  analyses  are  less  reliable  near  fronts  and  in  highly 
variable  areas. 

The  main  advantages  of  ICAPS  are  that  (1)  it  can  perform  a  synoptic 
water-mass  analysis  using  observed  data,  (2)  it  is  automated  to  the 
extent  that  a  nonoceanographer  or  nonacoustician  can  get  good  results 
from  it,  and  (3)  it  is  operational  and  available  now. 

ODSI  is  a  proposed  model,  yet  to  be  built,  that  will  provide  a 
variety  of  oceanographic  data  including  monthly  maximum,  minimum,  and 
"most  probable"  profiles  of  temperature,  salinity,  computed  sound 
speed,  and  sigma-T  from  surface  to  bottom.  The  spatial  grid  will 
typically  be  1°  x  1°  or  less,  depending  on  the  local  complexity. 

Coverage  will  be  worldwide.  The  model  will  also  provide  information  on 
the  frequency  of  current  and  frontal  boundaries  within  squares,  strong 
currents,  eddies,  "shelf  phenomena"  (not  explained),  and  mixed  layer 
thickness. 

ODSI's  disadvantage  is  that  it  remains  largely  conceptual  and 
untested.  We  have  seen  no  outputs  to  date,  and  more  information  on 
the  methodology  will  be  needed  for  an  accurate  evaluation.  The  developer 
estimates  that  the  analysis  for  the  entire  world  could  be  completed  in 
18  months  at  a  cost  of  ~  $180,000;  smaller  portions  could  be  provided 
on  a  prorated  basis.  We  are  unsure  of  the  accuracy  of  these  estimates. 

The  main  advantage  of  ODSI  is  the  variety  and  abundance  of  data 
it  promises  to  provide,  although  some  data  may  not  be  direct  J.y  appro¬ 
priate  to  LRAPP's  needs.  It  is  the  only  model  we  have  found  that 
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at tempts  to  deal  with  ocean  variability.  The  variable  spatial  grid 
will  be  an  advantage  in  areas  containing  fronts  and  currents.  The 
procedure  for  producing  the  profiles  will  atteupt  to  preserve  horizontal 
gradients  and  vertical  curvature. 

!!  I  MAX  is  designed  to  provide  deep  so^nd-speed  profiles  fe366  m) 
to  extend  an  observed  XBT-derived  sound-sp^ed  profile  to  the  bottom. 

The  basic  data  files  contain  annual  sound- speed  profiles  for  70  so-called 
homogeneous  areas  in  the  North  Atlantic,  69  areas  in  the  North  Pacific, 
md  50  areas  in  the  Indian  Ocean.  Spatial  resolution  for  data  selection 
is  presumed  to  be  1*  *  i°.  In  addition,  monthly  best-guess  sound-speed 
profiles  are  available  for  the  upper  366  m.  The  program  and  data  are 
computerized  for  shipboard  operation.  Profiles  are  a  Is,?  available  in 
atlas  reference  charts  by  area. 

SIMAS  is  unsuitable  because  its  use  of  homogeneous  area  profiles 
does  not  yield  realistic  horizontal  gradiencs,  so  it  cannot  produce 
reliable  inputs  for  range-dependent  models.  In  its  favor,  the  analysis 
is  complete  and  published  in  atlas  form.  Its  availability,  acquisition 
cost,  and  time  to  implement  are  unknown. 


Pi t  eti tiall  v  Sui  t ab  le  C and  t dates 

Is  a  CDC  6000  series  retrieval  system  for  selecting 
eti vl  r oilmen t a  1  ,..ir  meters  suitable  fui  input  to  NORDA’s  numerical 
acoustic  models.  The  retrieval  data  include  (1)  the  most  representa¬ 
tive  observed  seasonal  sound-speed  profile  (originally  produced  by 
RiiVP)  ,  by  ')*  *  V*  squares,  (2)  bathymetry  data  from  Seripps  Institution 
ot  Oceanography,  by  l*  *  1°  squares  or  tracks,  (3)  seasonal  significant 
wave  heights  by  5®  *  5*  squares,  and  (4)  FNVC  and  NAVOCEANO  bottom  loss 
classes  bv  l*  *  1*  squares.  Coverage  is  for  the  northern  hemisphere . 
Outputs  arc  on  great-circle  tracks.  Sound-speed  protiles  are  extended 
to  the  bottom  using  the  deep  profile  data  contained  in  KSVP ;  those 
data  are  the  average  sound  speeds  in  each  5*  *  5*  squat-  at  1000  a 
intervals  between  1000  a  and  10,000  a. 

The  main  advantage  of  AUTO-OCEAN  is  that  its  products  describe 
parameters  ttot  included  in  anv  other  model  (items  2,  3,  and  4  above); 
those  parameters  are  required  as  Inputs  ior  numerical  acoustic  models. 


i.-,  *■ 
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AUTO-OCEAN  also  operates  rapidly  and  inexpensively  (e.g.,  a  single 
season  of  parameters  for  a  50  nm  track  could  be  produced  in  ~  4  sec 
at  a  cost  of  ~  24  cents,  government  rate).  LRAPP  has  paid  for  the 
model,  and  it  is  available  now. 

AUTO-OCEAN's  5“  *  5*  spatial  resolution  for  the  sound-speed  model 
is  too  coarse,  and  the  model  lacks  temperature  and  salinity  data;  how¬ 
ever,  representative  sound-speed,  temperature,  and  salinity  profiles 
can  easily  be  Inserted  on  a  1*  *  1*  grid  spacing. 

GDEM  is  an  objective-analysis  program  that  produces  fields  of 
analyzed  temperature,  salinity,  and  sound  speed  on  a  30'  *  30'  grid 

It 

over  various  depth  ranges.  Preliminary  analysis  has  been  completed 
for  the  North  Pacific  and  Mediterranean  and  is  beginning  for  the  North 
Atlantic.  The  Indian  Ocean  segment  is  to  be  completed  sometime  in  the 
future.  Observational  data  are  insufficient  to  run  GDEM  in  the  southern 
hemisphere. 

The  (IDEM  component  for  the  upper  layer  models  seasonal  temperature 
between  the  surface  and  the  merge  depth  (~  400  m) annual  salinity 
between  the  surface  and  either  400  m  in  the  Mediterranean  or  800  a 
in  the  Pacific,  and  seasonal  sound  speed  from  the  surface  to  the 
merge  depth  in  both  oceans.  The  sound  speed  is  derived  from  the 
analyzed  temperature  and  salinity  fields. 

Below  the  surface  models  of  temperature ,  salinity,  and  sound  speed 
there  is  a  aiddepth,  two-season  (w I nier -summer)  sound-speed  model  from 
200  m  to  2450  ts.  It  is  being  augmented  in  the  North  Atlantic  by  a  mid¬ 
depth  temperature  and  salinity  model  (see  Appendix  C,  Attachment  3). 

The  lowest  component  is  an  annual  sound-speed  model  extending 
from  2000  »  to  the  bottom.  There  la  no  evidence  that  the  developer 
intends  to  replace  this  model  with  a  temperature  or  salinity  model  in 
the  near  future. 

This  taudvl  has  received  our  most  detailed  scrutiny  to  date. 

* 

Merge  depth,  a  variable  peculiar  to  each  ocean  basin,  is  the 
depth  at  which  two  adjacent  layers  ot  the  as.-wl.-i  arc  merged.  In  the 
earlv  19?V  version  of  the  model  for  the  North  Pacific,  for  example, 
that  depth  was  tixed  for  the  entire  basin  a*  400  m.  Ue  understand 
that  the  latest  versions  of  the  model  adjust  merge  depth  over  the 
basin  to  fit  the  observed  data. 
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Each  GDEM  component  uses  a  different  curve-fitting  technique  to 
obtain  the  best  fit  to  the  edited  observed  data.  Bathymetry  data  for 
tin*  North  Pacific  and  Mediterranean  are  on  a  30'  *  30'  grid. 

Of  all  Standard  Ocean  candidate*,  using  objective-analysis  tech¬ 
niques,  GDEM  comes  closest  to  reproducing  the  significant  features  in 
vertical  sound-speed  profiles.  Comparison  of  model  profiles  with 
typical  observed  profiles  for  the  North  Pacific  showed  that  GDEM  main¬ 
tains  horizontal  continuity  yet  preserves  frontal  structure  in  some  cases. 
It  appears  to  work  well  in  areas  with  sparse  data  except  the  southern 
hemisphere.  The  30;  *  30'  grid  spacing  is  the  smallest  of  any  model 
yet  evaluated.  Operation  is  rapid  (~  3  min  to  edit  an  entire  tape  for 
one  season  In  the  North  Pacific)  and  reported  to  be  simple. 

Several  problems  need  to  be  resolved  before  GDEM  could  become  the 
Standard  Ocean  retrieval  system.  One  is  a  tendency  toward  error  in 
modeling  the  upper  layers  at  the  higher  latitudes  in  the  i.orth  Pacific. 
The  developer  has  undertaken  a  redest,  of  that  portion  of  the  model 
and  may  already  have  solved  the  problem  (see  Appendix  C,  Attachment  3). 

A  second  problem  Is  an  inexplicable  seasonal  shift  in  the  GDEM  sound- 
speed  profiles  in  the  North  Pacific  near  30*N,  140*W,  between  300  and 
1600  m.  We  do  not  know  whether  the  recent  redesign  will  correct  that 
discrepancy.  Its  acoustic  significance,  however,  is  expected  to  be 
minimal.  The  redesigned  model  will  be  tested  in  PY  1980.  The  developer 
expects  GDEM  to  he  tully  operational  by  early  FY  1931. 


III.  CONCLUSIONS 


Of  the  eight  original  candidates  for  Standard  Ocean,  we  judge 
that  AUTO-OCEAN  and  GDEM  show  the  most  potential  for  meeting  LRAPP 's 
requirements,  and  in  a  timely  and  cost-effective  manner.  Both  models 
are  continually  being  evaluated  and  modified,  though  AUTO-OCEAN  is 
available  now.  Because  of  GDEM's  realistic  and  detailed  representa¬ 
tion  of  the  oceanic  vertical  and  horizontal  structure,  it  is  our  first 
choice  for  Standard  Ocean.  If  GDEM  is  not  available  for  installation 
by  early  FY  1981,  we  recommend  that  AUTO-OCEAN  be  installed  in  the 
interim  and  updated  using  typical  observed  profiles,  as  mentioned  in 
Sec.  II.  By  mid-FY  1980,  we  will  know  whether  GDEM  will  be  available 
on  schedule;  then  we  can  recommend  an  appropriate  course  of  action  to 
bring  a  Standard  Ocean  on-line  by  early  FY  1981.  When  GDEM  is  approved 
and  released,  we  recommend  that  its  output  be  stored  in  the  AUTO-OCEAN 
framework.  The  advantage  is  that  AUTO-OCEAN  is  a  data  retrieval  system 
known  to  work;  it  should  not  be  too  difficult  to  replace  AUTO-OCEAN's 
outmoded  sound -speed  profiles  with  the  new  GDEM  profiles.  Storing 
GDEM  output  in  the  AUTO-OCEAN  framework  would  enable  the  acoustic 
modeler  to  retrieve  sound  speed,  bottom  loss  class,  bathymetry,  and 
wave  height  in  a  single  model  run  instead  of  the  two  or  three  runs 
otherwise  required. 

The  only  LRAPP  requirement  no  candidate  meets  is  the  one  for  ocean 
variability .  We  suggest  that  LRAPP  sponsor  a  study  of  ocean  variability 
to  produce  a  detailed  1°  *  1°  variability  index  to  be  included  in 
Standard  Ocear..  We  do  not  expect  that  the  index  could  be  implemented 
in  t Lme  for  the  initial  installation  of  Standard  Ocean  in  early  FY  1981. 


i 
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Appendix  A 


AUTO-OCEAN 


AUTO-OCEAN  is  a  CDC  6000  series  retrieval  system  of  environmental 
parameters  suitable  for  input  in  numerical  acoustic  models.  It  is 


operational  at  Eglin  AFB  and  is  accessed  via  a  remote  terminal  at  NORDA 
(Code  320) .  AUTO-OCEAN  provides  the  following  environmental  informa¬ 
tion  for  the  entire  northern  hemisphere:  (1)  bathymetry  (Scripps 
Institution  of  Oceanography)  per  1°  square  or  track,  (2)  seasonal  sig¬ 
nificant  wave  heights  per  5°  square,  (3)  bottom  loss  classes  (FNWC  1-5 
and  N00  1-9)  per  1°  square,  and  (4)  seasonal  vertical  sound-speed  pro¬ 
files  (using  RSVP)  per  5°  square.  The  sound-speed  profiles  are  the 
most  representative  observed  for  the  given  square  and  season.  The 
profiles  are  extended  to  the  bottom  using  RSVP's  deep  profile  data, 


which  are  average  sound-speed  values  for  each  5°  <  5°  square  at  1000  m 


intervals  between  1000  m  and  10,000  m. 

The  environmental  information  is  generated  along  a  great-circle 


path  given  either  an  initial  point  (lat,  Ion),  bearing,  and  maximum 
range,  or  two  points  (initial  and  final  lat.  Ion).  The  parameters  are 
retrieved  at  each  point  where  the  selected  track  intersects  a  1*  lat. 
Ion  grid  on  the  earth's  surface.  The  data  sources  are  the  NODC  ocean 
station  data  file  to  1972,  FNWC  and  NAVOCEANO  bottom  loss  classes,  AES0 
interim  wave  height  data  bank,  and  the  Scripps  bathymetric  data  file. 

Turnaround  tune  for  AUTO-OCEAN  is  typically  4  sec  for  a  50  tin 
track;  cost  is  about  24  cents  for  one  season.  The  system  is  available 
now.  Us  technical  quality  is  acceptable  for  most  acoustic  model  runs. 


*Thls  appendix  draws  on  Information  provided  in  a  letter  to  Daub in 
from  Hashlmoto,  28  March  1979. 

See  J.  J.  Audet,  Jr.,  and  0.  E .  Vega,  AaN."  rlv/:.*- 

r<t- tripl'd!  t,*"  ('b'.'TT 1 ,  Acoustic  Environmental  Support  Detachment, 

ONR,  AESD  Technical  Note  TN- 74-03,  October  1974;  and  Audet,  Sauaowi 
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depending  on  application  and  purpose.  LRAPP  has  already  paid  for  the 
product,  so  It  could  be  acquired  at  no  cost. 

If  AUTO-OCEAN  is  to  serve  as  an  interim  or  permanent  Standard 
Ocean,  we  recommend  that  the  current  5°-square  sound-speed  profiles 
be  replaced  by  representative  l°-square  sound-speed,  temperature,  and 
salinity  profiles  that  are  compatible  with  the  1°  bathymetry  and  1° 

bottom-loss-class  data  files.  Those  profiles  could  be  produced  by 

★ 

selecting  the  most  typical  observed  profiles  or  by  an  objective  analy¬ 
sis  technique  (e.g.,  GDEM — see  Appendix  C) .  Objective  analysis  might 
preserve  the  horizontal  continuity  of  the  field  better  than  observed 
prof  lies. 

* 

For  example,  J.  G.  Colborn  and  J.  D.  Pugh,  A  Pi'oceduve  for  Selec¬ 
tion  rf  Typical  Sound  Speed  Pro f ilea ,  Naval  Undersea  Center,  NUC  TN 
1006,  May  1973. 
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Appendix  B 
FIB/EOTS/EXTRA* 

Hie  FIB/EOTS  (Field  by  Information  Blending/Extended  Ocean  Thermal 
Structure)  retrieval  system  is  used  at  Fleet  Numerical  Weather  Central 
(FNWC)  for  predicting  oceanic  thermal  structure.  The  system  now  oper¬ 
ates  on  a  CDC  6500  computer;  FNWC  will  use  a  Cyber  175  in  the  future. 
FIB/EOTS/EXTRA,  discussed  below,  is  the  FIB/EOTS  version  that  produces 
sound-speed  profiles. 

ANALYSIS 

FIB/EOTS  accepts  irregularly  spaced  (temporal  and  spatial)  tem¬ 
perature  profiles  as  input.  Those  single  observations  are  subjectively 
weighted  according  to  relative  accuracy  and  interpolated  to  the  nodes 
of  FNWC 1 s  standard  63  *  63  grid  of  the  northern  hemisphere  (see  Fig. 
B.l).  At  60°N ,  that  grid  has  a  spatial  resolution  of  Ax  ~  Ay  ~  380  km. 
For  limited  geographical  areas,  a  higher  resolution  grid  is  available 
down  to  1/8  the  standard.  Then,  "fields,"  horizontally  analyzed  sur¬ 
faces  of  temperature  and  its  first  and  second  derivatives,  are  fitted 
to  the  data  values  at  the  grid  nodes.  The  fields  are  computed  at  26 
depths  (including  standard  depths)  between  the  surface  and  5000  m,  with 
Increased  resolution  near  the  primary  layer  depth  (i.e.,  the  point  of 
the  largest  vertical  curvature  in  the  temperature  profile). 

for  regions  where  no  observed  data  are  available,  FIB/EOTS  sub¬ 
stitutes  archived  climatology  values,  though  not  limited  to  the  XBT 
(T-4)  data  format,  FIB/EOTS  uses  primarily  XBT  observations  as  Input. 

On  a  typical  day,  400  are  received,  most  as  fleet  messages  (20-30  tem¬ 
perature-depth  pairs)  rather  than  actual  XBT  traces.  For  an  input  of 
that  size,  the  turnaround  time  is  approximately  1/2  hr.^ 

* 

This  appendix  draws  on  information  provided  in  a  letter  to  Daubin 
from  Ryan.  2  April  1479. 

♦ 

FIB  objective  methodology ,  the  basis  of  FIB/EOTS,  is  further  de¬ 
scribed  in  the  attachment  to  this  appendix. 


tandard  63  x  63  grid  of  the  northern  hemisphere 
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OUTPUT  PRODUCTS 

Temperature  Profiles 

Temperature  profiles  are  computed  at  each  grid  node  by  merging 
the  FIB-analyzed  dynamic  upper  levels  (£400  m)  with  the  3tatic  deeper 
level  (>400  m)  climatology  in  a  way  to  minimize  an  error  functional. 

The  error  functional  is  essentially  a  least-square  minimization  of 
the  gradient  and  curvature  of  the  computed  profile.  The  merging  pro¬ 
cess  tends  to  shift  the  upper  portion  of  the  deeper  level  climatology 
toward  the  XBT  observations  near  400  m.  Any  climatological  anomalies 
in  the  dynamic  upper  levels  will  be  reflected  in  the  profiles  at  greater 
depths.  Because  of  the  statistical  nature  of  the  FIB  analysis,  hori¬ 
zontal  gradients  will  diminish  unless  the  spatial  data  have  been  sampled 
adequately.  However,  the  user  can  impose  boundary  conditions  on  the 
FIB  analysis  to  retain  certain  gradient  features. 

Sound-Speed  Profiles 

The  FNWC  program  EXTRA  (EXTRAction)  must  be  added  to  FIB/EOTS 
to  produce  sound-speed  profiles.  Using  the  FIB/EOTS  temperature  fields 

and  unanalyzed,  archived  salinity  data,  EXTRA  computes  sound  speed 

* 

from  Leroy  s  equations.  The  program  does  not  check  the  vertical  sta¬ 
bility  of  the  computed  temperature/salinity  values  at  the  grid  nodes. 
Sound-speed  profiles  at  points  other  than  grid  nodes  are  computed  by 
linear  interpolation  between  adjacent  nodal  values. 

ADVANTAGES 

FiB/EOTS  has  a  number  of  positive  qualities  for  oceanic  tempera¬ 
ture  analysis: 

•  It  is  an  operational  program. 

•  Its  documentation  and  source  code  are  available. 

•  It  is  continuously  being  tested  and  evaluated. 

• 

C.  C.  Leroy,  "Development  of  Simple  Equations  for  Accurate  and 
More  Realistic  Calculation  of  the  Speed  of  Sound  in  Seawater,"  J. 

Aocuet.  Soc.  Am.,  Vol.  46,  1969,  pp.  216-226. 
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•  It  is  highly  automated. 

a  The  methodology  is  internally  consistent. 

•  It  has  the  capability  of  variable  temporal  resolution. 

•  It  easily  accepts  new  observations. 

•  Error  estimates  are  made  on  the  finai  output  fields. 

•  It  adequately  preserves  layer  depth. 

•  It  does  not  overburden  computer  resources. 

•  The  FI B/EOTS /EXTRA  objective  analysis  technique  may  be  suit¬ 
able  for  generating  some  Standard  Ocean  data  sets. 

DISADVANTAGES 

For  adaptation  as  Standard  Ocean,  FI B/EOTS /EXTRA  has  several 
weaknesses : 

•  The  spatial  resolution  of  the  standard  FNWC  63  *  63  grid  is 
too  coarse. 

•  Only  temperature  is  analyzed;  the  salinity  values  are  archived 
ones . 

•  The  analysis  is  sensitive  to  the  scarcity  of  synoptic  XBT  ob¬ 
servations  and  thus  relies  heavily  on  static  climatology. 

•  Vertical  profile  integrity  is  compromised  by  the  merging  of 
shallow  (<400  m)  synoptic  XBT  observations  and  deeper  clim¬ 
atology. 

•  The  climatology  fields  (temperature  and  salinity)  used  as 
inputs  to  the  FIB  analysis  are  of  questionable  suitability. 
(This  weakness  may  be  corrected  when  FNWC  replaces  the  cur¬ 
rent  climatology  with  the  new  and  larger  MOODS  data  set,  as 
planned. ) 
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Attachment 

FIB  OBJECTIVE  ANALYSIS  METHODOLOGY* 

The  Field  of  Information  Blending  (FIB)  is  a  powerful  analysis  tech¬ 
nique  applicable  to  virtually  any  two-dimensional  physical  variable.  As 
of  1973,  FIB  programs  were  being  used  at  FNWC  to  analyze  sea  surface 
temperature  and  sea  level  pressure  distributions.  As  of  June  1975,  an 
FIB  adaptation  for  wind  analysis  (u,  v)  has  been  formulated  and  other 
applications  are  under  development. 

The  FIB  technique  analyzes  the  distribution  of  a  variable  by  blend¬ 
ing  measurements  of  the  variable  and  its  gradients,  which  come  from 
different  sources  and  locations.  The  program  uses  reports  from  various 
observation  stations,  with  estimates  of  reliability,  and  it  accepts 
regional  or  whole  field  estimates  of  the  parameter  and  its  derivatives 
(gradient,  Laplaeian,  etc.).  It  checks  all  input  data,  rejects  gross 
errors,  and  assembles  the  data.  From  this,  it  blends  or  analyzes  to 
produce  the  optimum  analysis  which  best  fits  all  the  information  at  hand. 
The  technique  also  produces  grid-point  reliabilities  of  the  final  pro¬ 
duct.  All  input  data  are  reevaluated  individually  by  comparison  with 
tiie  blended  analysis,  which  includes  the  interacting  effects  of  all 
information  that  went  Into  the  analysis. 

Reliability  or  weight  is  a  measure  of  the  worth  of  a  piece  of  in¬ 
formation.  in  every  step  of  the  FIB  process,  information  exercises  only 
the  degree  of  influence  warranted  by  its  reliability  value  at  that  par¬ 
ticular  stage  of  the  analysis.  The  ability  to  compute  the  information's 
reliability  is  the  key  to  the  FIB  technique. 

FIB  has  six  component  operations,  to  be  discussed  in  the  following 
sections . 


This  attachment  is  paraphrased  from  U. S.  Naval  Weather  Service  Nu¬ 
merical  Environmental  h'oducts  Manual,  NAVAIR  50-IG-522,  1  June  1975,  pp. 
3.10-1  and  3.10-2.  According  to  Jack  Kaitala  (FNWC),  the  description 
applies  to  the  current  FIB  analysis. 
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FIRST-GUESS  FIELD  PREPARATION  OF  INITIALIZATION 

Thi'  first  guess  Is  an  estimate  of  what  the  analysis  will  be  with¬ 
out  considering  current  data.  It  provides  continuity  into  data-sparse 
areas  and  gives  an  estimate  of  the  shape  (gradients,  curvature,  etc.) 
of  the  field.  In  sparse-data  areas,  the  accuracy  of  the  first  analysis 
depends  partly  upon  first-guess  accuracy.  The  first  guess  is  also  useful 
for  keeping  "impossible"  data  from  entering  the  analysis  by  indicating 
the  approximate  values  expected  in  an  area.  Information  is  thus  tested 
for  credibility  against  the  first  guess. 

The  first  guess  in  objective  analysis  models  is  either  a  previous 
analysis,  extrapolated  to  analysis  time,  or  a  prognostic  chart  verified 
at  analysis  time. 

FIB  has  the  unique  capability  of  accepting  several  first-guess 
fields,  each  weighted  by  its  proportionate  value.  Later  in  the  program, 
KLB  can  individually  reevaluate  the  worth  of  each  first-guess  field. 

If  .1  previous  analysis  is  to  be  one  of  the  first-guess  fields,  the  FIB 
program  has  a  special  steering  subroutine  to  bring  the  old  analysis  up 
to  analysis  time.  In  the  sea-level  pressure  application,  the  500-mb  SR 
(residual)  height  field  steers  the  previous  analyses  the  appropriate 
distance  and  direction.  This  process,  called  kineraatical  extrapolat ion , 
gives  a  conservative  first  guess  and  is  very  informative  when  used  with 
a  more  sophisticated,  primitive-equation  forecast. 

ASSEMBLY  OF  NEW  INFORMATION 

Reports  of  the  parameter  being  analyzed  are  placed  at  their  proper 
geographical  positions.  The  first-guess  field  value  is  Interpolated  at 
the  report  location,  and  the  values  are  compared.  In  sene  analysis  sys¬ 
tems,  only  an  arbitrarily  assigned  difference  is  allowed.  If  the  allowed 
difference  is  exceeded,  the  report  is  thrown  out.  In  the  FIB  method, 
the  difference  varies  with  the  magnitude  of  the  gradients  near  the  report 
and,  for  some  parameters,  with  latitude. 

If  the  report  passes  the  gross-error  check,  it  is  assigned  a  re¬ 
liability  or  weight.  This  is  based  on  the  standard  deviation  of  the 
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errors  associated  with  that  type  of  report — the  larger  the  standard 
deviation,  the  smaller  the  reliability.  Other  factors  involved  iri  deter¬ 
mining  the  reliability  vary  with  the  parameter;  examples  include  magni¬ 
tude  of  the  gradients  near  the  report,  age  of  the  report,  and  station 
elevation . 

The  difference  between  the  interpolated  first-guess  value  and  the 
report  value  is  applied  at  the  nearest  grid  point  as  a  correction  to  the 
grid-point  guess  value.  The  total  assembled  value  of  the  parameter  at 
a  grid  point  is  a  weighted  mean  of  all  the  data  referred  to  the  grid 
point,  with  each  value  contributing  to  the  mean  in  accordance  with  its 
re  1 iabi 1 i ty . 


BLENDING  FOR  THE  PARAMETER 

Blending  is  the  analysis  stage,  corresponding  to  the  drawing  of 
isolines  by  a  hand  analyst.  The  assembly  step  combined  reports  at  their 
nearest  grid  points.  The  grid-point  information  is  now  spread  to  sur¬ 
rounding  grid  points  by  reference  to  previously  derived  gradients  and 
higher  order  fields.  The  degree  of  spreading  is  increased  with  higher 
reliability  in  the  gradient  and  other  spreading  fields. 

After  blending,  each  grid  point  will  have  a  new  parameter  value, 
reflecting  surrounding  information  as  well  as  information  at  the  grid 
point  itself.  This  is  an  optimum  compromise  of  all  the  weighted  infor¬ 
mat  ion . 


RELIABILITY  HELD  OF  THE  BLENDED  PARAMETER 

Next,  a  reliability  field  Is  computed  for  the  blended  parameter 
in  preparation  for  the  next  step  of  reevaluation  and  error  checking. 

The  blending  process  results  in  new  parameter  values  at  each  grid 
point,  and  these  values  will  have  a  different  reliability  or  weight. 
For  example,  if  one  grid  point  had  nothing  but  first-guess  information 
before  blending,  its  reliability  would  Le  much  lower  than  surrounding 
grid  points  that  have  information  from  -.everal  observations.  Blending 
spreads  the  information  from  the  high-reliability  grid  points  to  the 
ones  chat  have  lower  reliability.  Blending  increases  reliability  at 
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n  I  I  grid  points,  reflecting  the  additional  information  flowing  in  from 
surrounding  grid  points.  Even  the  low-reliability  points  can  add  infor¬ 
mation  to  surrounding  areas.  The  interaction  between  the  grid  points 
is  greatest  over  a  one-grid  interval  and  diminishes  with  distance.  The 
strength  of  interaction  is  limited  by  the  gradient  weight,  i.e.,  if  the 
gradient  is  known  only  slightly  (low  weight),  even  adjacent  grid  points 
will  have  little  effect  on  each  other. 

KEEVA1.UATI0N  AND  LATERAL  REJECTION 

FIR  uses  the  blended  parameter  field  and  grid-point  reliabilities 
to  reevaluate  each  piece  of  information  that  entered  into  the  analysis. 
The  reevaluation  provides  quality  measures  for  each  observation  and 
each  first-guess  field.  The  analysis  cycle  will  be  repeated  using  the 
reweighted  information.  TVie  reevaluacion  stage  is  a  vital  and  integral 
part  of  the  FIB  technique.  This  allows  a  second  or  even  third  analysis 
pass  with  ever- improving  weights. 

To  reevaluate  reports  of  any  parameter,  a  statistical  measure  is 
computed  for  each  report.  This  measure  indicates  the  accuracy  of  a  re¬ 
port  compared  with  the  accuracy  expected  from  the  designed  reliability. 
Each  report,  with  its  weight,  is  individually  removed  from  its  grid 
point  and  compared  with  what  remains,  or  the  '‘background."  If  the  re¬ 
port  is  within  its  expected  error,  no  change  is  made  in  its  reliability, 
if  the  error  is  greater  than  expected  but  within  some  upper  limit,  the 
report’s  reliability  is  reduced.  if  the  error  limit  5s  exceeded,  the 
report  is  rejected  (l.e.,  its  weight  becomes  zero)  and  It  will  have  no 
effect  whe:.  the  next  assembly  and  blending  is  made. 

The  reliability  of  the  first-guess  fields  can  be  similarly  evalu¬ 
ated.  If  new  ini  ,'rtaat ion  disagrees  with  the  first  guess,  the  weight 
of  the  latter  is  revised.  In  some  app 1 icat lone  involving  rapid  change, 
such  as  sea-level  pressure,  the  first-guess  weight  is  so  small  that 
reevaluation  is  not  necessary. 

FEASALYStS 

The  rean4lysls  begins  by  returning  •  >  the  asscnbly  stage.  Hie  new 


-29- 


assembly  starts  with  the  first-guess  fields,  which  may  be  reweighted, 
and  the  reports  are  assembled  with  their  reevaluated  weights.  The  whole 
cycle  is  repeated  exactly  as  before.  In  the  final  pass  (second  or  third 
analysis),  the  program  skips  the  reevaluation  and  proceeds  to  the  output 
section.  The  final  analysis  is  stored  in  the  computer  for  transmission 
and  for  input  to  other  programs. 
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Appendix  C 
GDEM* 


The  Generalized  Digital  Environmental  Model  (GDEM)  being  developed 

by  Thomas  Davis  of  NAVOCEANO  is  an  objective  analysis  program  designed 

to  produce  seasonal  fields  of  analyzed  temperature,  salinity,  and  sound 

speed  on  a  30'  *  30'  grid.  The  model’s  output  can  be  adapted  for  use 

in  an  automated  environmental  data  retrieval  program  to  provide  range- 

dependent  inputs  for  acoustic  modeling.  GDEM  is  not  yet  operational, 

and  the  modeling  processes  are  continually  being  evaluated  and  re- 

-f- 

developed.  Therefore,  little  documentation  is  available,  and  current 
descriptions  must  be  regarded  as  tentative. 

PRELIMINARY  EVALUATION 

To  evaluate  GDEM's  ability  to  reproduce  the  typical  vertical 
structure  in  selected  oceanographic  regions,  NOSC  selected  real  sound- 
speed  profiles  to  represent  six  Northeast  Pacific  sound-speed  provinces 
for  the  summer  season  and  compared  them  with  model  outputs.  We  re¬ 
ported  the  results  to  LRAPP  on  30  April  1979  (reproduced  as  Attachment 
2).  Some  province  profiles  matched  well,  confirming  GDEM's  potential 
ability  to  reproduce  real  oceanographic  vertical  structure.  , 

The  results  also  uncovered  a  technical  problem  regarding  the  model¬ 
ing  of  shallow-channel  sound-speed  structures  observed  around  latitude 
49 °N.  In  some  provinces  the  difference  between  the  "typical"  rj>.al  arid 
GDEM  profiles  exceeded  the  standard  deviation  for  all  real  observations 
in  the  province  at  a  number  of  depths.  The  greatest  mismatches  occurred 


This  appendix  was  written  by  Colborn,  who  also  analyzed  the  NOSC 
comparison  results, 
j* 

NAVOCEANO  described  the  basic  techniques  of  an  early  version  of 
GDEM  in  an  informal,  unpublished  report  (reproduced  as  Attachment  1  to 
this  appendiv'i  .  A  mote  detailed  mathematical  interpretation  of  the 
model,  "LRAPP  Objective.  Analysis,"  Ocean  Data  Systems,  Inc.,  20  March 
1979,  unpublished  report,  has  been  outdated  by  later  model  revisions. 
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;it  230-300  m.  Tn  one  instance,  the  model  profile  shifted  the  main 
channel  ax  L.s  to  a  depth  outside  the  range  of  observed  values  near  the 
test  location. 

Those  discrepancies  prompted  NAVOCEANO  to  redesign  portions  of 
the  model.  At  first  the  problem  was  thought  to  lie  in  the  merging  of 
the  upper  Butterworth  filter  model  and  the  middepth  orthogonal  model. 
NAVOCEANO  later  found  that  replacing  the  quadratic  tail  in  the  upper- 
layer  temperature  model  by  an  exponential  function  eliminated  the 
shallow-channel  discrepancies  (see  Attachment  3  for  a  recent  memorandum 
reporting  on  the  GDEM  redesign) . 

RECENT  EVALUATIONS 

Vertical  Structure 

While  awaiting  the  results  of  the  model  redesign,  NOSC  continued 
comparing  data  from  the  original  model  with  North  Pacific  observations. 
'Ihe  earlier  comparisons  used  real  profiles  from  an  analysis  defining 
sound-speed  provinces.  Those  profiles  were  developed  to  represent  a 
large  region,  whereas  GDEM  is  designed  to  produce  a  profile  representa¬ 
tive  of  a  single  location.  To  match  that  design  characteristic,  we 
decided  in  the  later  comparisons  to  select  real  profiles  representing 
point  locations.  Because  a  reasonable  number  of  observations  are 
required  for  selecting  a  typical  profile,  the  data  inventory  W3S 
searched  to  locate  1°  squares  with  high  data  density  for  evaluation. 
Figure  C.l  shows  the  location  of  the  six  1°  squares  selected.  Locations 
1,  2,  3,  and  5  met  the  criterion  of  a  reasonable  number  of  observations 
for  winter  ( January-March)  and  summer  (July-September) .  Location  4  met 
the  criterion  only  for  the  summer,  and  location  6  did  not  provide  ade¬ 
quate  data  for  either  season. 

Figures  C.2  through  C.13  graphically  depict  the  comparisons  between 
the  selected  typical  profile  for  the  specified  1°  square  and  the  GDEM 
profile  for  the  central  coordinates  of  the  square.  The  maximum  and 
minimum  profile  envelopes  for  all  observed  data  are  also  plotted. 

Figures  C.14  through  C.25  overplot  the  GDEM  iprofile  on  the  composite 
plot  for  all  deep  observed  profiles  in  each  seasonal  l°-square  data 
set . 


Locations  of  6DEM- typical  observed  profile  comparison 
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Fig.  C.6--GDEM-typical  comparison,  location  3  (34°N,  164°E),  winter 
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Fig.  C.9--GDEM-typical  comparison,  location  4  (51*N,  156*E),  summer 
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Location  1.  Figures  C.2  and  C,3  indicate  that  the  shapes  of  the 
typical  profile  and  the  GDEM  output  are  similar  in  winter  and  summer. 
However,  comparison  of  just  the  GDEM  profiles  for  winter  and  summer 
indicates  a  large  sound-speed  difference  extending  to  a  depth  of 
2000  m.  That  seasonal  difference  is  not  found  in  the  typical  profiles 
or  in  the  ouserved  data  sets  containing  more  than  400  observations 
for  each  season.  Figures  0.14  and  C.15  also  indicate  that  the  GDEM 
middepth  sound  speeds  fall  positively  outside  the  data  envelope  during 
the  winter  and  negatively  outside  during  the  summer.  That  variation 
in  the  GDEM  output  cannot  be  readily  explained  by  known  oceanographic 
seasonal  variability  for  the  area.  The  potential  problem  appears  to 
lie  in  the  middepth  sound-speed  portion  of  the  model;  further  investiga¬ 
tion  is  needed  to  make  certain. 

Location  2.  The  structural  comparison  for  winter  (Fig.  C.4)  looks 
reasonable,  though  the  GDEM  profile  falls  outside  the  observed  data 
envelope  on  the  positive  side.  This  is  also  apparent  in  the  overplot 
of  the  GDEM  profile  on  the  composite  profile  (Fig.  C.16).  The  summer 
profile  (Fig.  G.5)  appears  similar  to  the  typical,  but  a  portion 
deviates  positively  from  the  observed  data  set  (also  indicated  in  Fig. 

C .  1 7) .  These  profile  mismatches  resemble  the  mismatches  encountered 
in  the  earlier  comparisons,  so  the  redesigned  upper-layer  model  may 
correct  both.  We  expect  to  evaluate  the  new  GDEM  when  model  output  is 
available  in  FY  1980. 

Loc at  Ion  1.  Figures  C.b  and  C.7  indicate  a  reasonably  good  match 
in  the  comparisons.  Both  CDEM  and  typical  profiles  indicate  seasonal 
deviations  at  depths  below  the  normal  seasonal  limit  of  100-200  m.  An 
explanation  is  suggested  by  the  composite  plot  for  the  summer  (Fig. 
C.19).  The  historical  transit  of  the  southern  part  of  the  eubtropical- 
subarctic  transition  zone  past  location  3  is  apparent  in  the  bimodal 
structure  ot  the  upper  thermocline.  The  GDEM  model  produces  an  inter¬ 
mediate  -.tructure  while  tiie  selection  procedure  for  the  typical  profile 
targets  the  shallow  thermocline  data  subset. 

Location  4.  The  selection  of  location  4  was  an  attempt  to  check 
the  CDEM  model  In  seml-isoiated  shallow  water.  However,  insufficient 
data  In  the  winter  (see  Fig.  C. 20)  and  high  variability  in  the  summer 
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(Fig.  C.21)  precluded  a  meaningful  model  comparison.  The  evidence  in 
Figs.  C.8  and  C.9  suggests  that  the  model  predicted  quite  well  the 
general  sound-speed  structure  for  this  location. 

Location  5.  Location  5  is  also  in  an  area  of  high  variability, 
as  the  composite  plots  in  Figs.  C.22  and  C.23  show.  The  GDEM  output 
profile  compares  quite  favorably  with  the  typical  and  composite  data 
during  the  summer.  The  winter  structure  is  similar,  but  Figs.  C.8 
and  C.22  show  a  slight  negative  shift  of  the  profile.  No  explanation 
is  apparent. 

Location  6.  Location  6  comparisons  are  inconclusive  because  of 
the  lack  of  observational  data.  The  winter  season  contained  only  one 
observation,  and  few  data  were  available  for  the  summer.  The  GDEM- 
composite  profile  overlays  in  Figs.  C.24  and  C.25  show  generally 
similar  structures.  A  tendency  toward  mismatch  is  seen  in  the  depth 
zone  where  the  model  is  being  redesigned.  The  new  model  will  be 
evaluated  for  location  6  in  FY  1980. 

In  summary,  GDEM  exhibits  the  same  potential  to  reasonably  model 
the  oceanic  sound-speed  structure  for  acoustic  modeling  that  it  did 
in  the  earlier  comparisons.  The  problem  remaining  with  upper-structure 
modeling  at  higher  latitudes  may  already  have  been  solved  by  the  re¬ 
design.  We  do  not  know  whether  the  redesign  will  correct  the  second 
problem  noted,  the  apparent  seasonal  shift  observed  in  the  CDEM  pro¬ 
files  for  location  i.  That  anomaly  does  not,  however,  seriously 
impair  the  model's  potential  for  reasonably  reproducing  vertical  struc¬ 
ture  and  acoustic  . ransmtsslon.  Further  evaluation  in  FY  1980  will 
reveal  more  about  the  seasonal -sh l f t  phenomenon. 

Horizontal  tni_e&r It v 

For  Standard  Ocean  application,  it  is  important  that  GDEM  pre¬ 
serve  the  horizontal  integrity  of  the  structure  for  range -dependent 
acoustic  modeling.  Judging  hv  the  development  procedures  used  to 
produce  GDEM,  preservat iou  of  horizontal  Integrity  will  be  a  funda¬ 
mental  asset  of  the  model. 

NAVOCEANO  gave  NOSC  a  limited  opportunity  to  test  this  charac¬ 


teristic  by  providing  a  single  GDEM-produced  sotmd-speed  contour 
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for  evaluation.  The  contour  pertained  to  a  section  of  the  great- 
circle  track  (Fig.  C.26)  extending  from  30°  to  45°N  through  the  Kuroshio 
and  Oyashio  frontal  zones  in  the  Northwest  Pacific  (Fig.  C.27).  Archi¬ 
val  observed  data  have  not  been  processed  for  comparison  with  the  GDEM 
section,  but  data  in  the  literature  allow  indirect  evaluation.  The 
key  aspects  of  model  output  to  consider  are  the  frontal  zone  gradients 
and  the  absolute  sound-speed  values. 

We  used  two  synoptic  sound-speed  sections  reported  by  Roden  for 
comparison.  The  first  is  a  north-south  section  at  168°E  and  at  the 
same  latitudes  as  the  GDEM  section,  for  April.  Although  located  con¬ 
siderably  east:  of  the  GDEM  section,  the  Roden  section  shows  a  horizontal 
gradient  of  24  m/s/l°  lat  at  200  m  depth  at  the  subarctic  front  that 
compares  with  the  GDEM  section  at  the  Oyashio  frontal  zone.  It  is  dif¬ 
ficult  to  distinguish  contours  and  to  compare  structure  in  the  upper 
200  m  on  the  GDEM  plot,  and  the  season  it  represents  is  unknown.  The 
second  Roden  section  at  154°E  for  April  provides  a  very  good  compari¬ 
son.  *  The  gradient  across  the  Kuroshio  front  in  this  Roden  section  is 
approximately  28  m/s/l°  lat  at  200  m  depth,  while  the  GDEM  gradient  is 
somewhat  less  at  15—20  m/s/1  lat.  Across  the  Oyashio  fiont  the  grad 
lents  are  very  similar  at  approximately  25  m/s/l°  lat. 

These  remarkable  results  indicate  that  smoothing  and  editing 
techniques  used  in  the  preprocessing  of  data  for  GDEM  do  not  obliterate 
significant  high-gradient  structural  features  such  as  large  frontal 
zones.  Tn  FY  1980,  more  comprehensive  range-dependent  testing  of  the 
model  should  be  performed  to  evaluate  the  horizontal  field  output. 

*G.  I.  Roden,  "Temperature  and  Salinity  Fronts  at  the  Boundaries 
of  the  Subarctic-Subtropical  Transition  Zone  in  the  Western  Pacific," 

J.  Geophys.  Res.}  Vol.  77,  No.  36,  1972,  pp.  7175-7187. 

V..  I.  Roden,  "On  North  Pacific  Temperature,  Salinity,  Sound 
Velocity  and  Density  Fronts  and  Their  Relation  to  the  Wind  and  Energy 
Flux  Fields,"  J.  Phys.  Oceano Vol.  5,  No.  4,  1975,  pp.  557-571. 


Fig.  C.26 — Great-circle 


contour  of  sound-speed  field,  30-45  N 
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Attachment  1 

DESCRIPTION  OF  TECHNIQUES  USED  FOR  DEVELOPMENT  OF 
SOUND  SPEED  PROVINCES  FROM  THE  NAVOCEANO  3D  MODEL 

INTRODUCTION 

The  data  set  used  in  deriving  the  Province  Chart  is  actually  a 
three-dimensional  seasonal  model  of  sound  speed.  The  model  consists 
of  three  parts:  a  surface  model  of  the  upper  400  meters  consisting 
of  annual  salinity  and  seasonal  temperature  (3  months);  a  two  season 
sound  speed  model  from  200m  to  2450m;  and  a  one  season  deep  sound 
speed  model  from  2000m  to  the  bottom.  A  bathymetry  file  is  also 
included . 

DATA  BASE 

The  basic  data  base  is  our  Ocean  Station  File  (OSTA)  which  con¬ 
tains  all  available  Nansen  cast  data.  This  file  was  used  for  all 
parameters.  To  better  define  the  near  surface  model  seasonally,  a 
XBT  file  supplemented  the  temperature  portion  of  the  model. 

DATA  ED  IT  INC 

There  are  two  phases  to  the  data  editing — automatic  and  manual. 

The  automatic  method  consists  of  a  computer  routine  which  does 
the  following: 

1.  Breaks  the  region  into  small  squares. 

2.  Computes  mean  and  standard  deviation  of  data  within  each 
square . 

3.  Eliminates  all  data  outside  i  one  standard  deviation  from  mean. 

In  practice,  a  double  pass  is  made  shifting  the  squares  to  better  eval¬ 
uate  stations  which  may  have  been  on  the  border.  For  deep  sound  speed 
data  a  5  degree  square  is  used  and  each  level  is  evaluated  separately. 

Unpublished  technical  description  by  NAVOCEANO,  Code  3300, 
spring  id 79. 
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For  data  which  consists  of  a  set  of  coefficients  for  each  station, 
a  different  approach  is  used.  The  profile  is  derived  from  the  coef¬ 
ficients  and  a  mean  profile  for  the  entire  data  set  is  computed.  All 
stations  are  compared  with  this  mean  profile  deriving  an  unnortnali zed 
coherency  coefficient  for  each  station.  These  coherency  coefficients 
are  then  used  as  the  values  in  editing.  A  two-degree  square  is  used. 

At  this  point  season  becomes  important  and  it  is  at  this  stage  that 
the  seasons  are  separated.  A  separate  run  is  therefore  required  for 
each  season. 

The  manual  procedure  consists  of  contouring  the  data,  locating 
questionable  areas,  listing  the  stations  within  and  near  these  areas, 
manually  selecting  the  bad  stations  and  running  a  cleanup  routine  to 
create  a  new  data  tape  without  the  bad  stations. 

CURVE  FITTING  METHODS 

The  heart  of  the  basic  model  is  the  orthogonal  polynomial  least 
squares  fit  to  each  input  profile  of  sound  speed.  One  of  the  problems 
with  tiie  orthogonal  polynomial  is  that  very  large  changes,  as  occur 
near  the  sea  surface,  are  not  only  poorly  fit,  but  also  cause  a  rippling 
effect  down  the  profile.  For  this  reason  the  upper  200m  were  not 
originally  used.  This  curve  fitting  technique  is  used  on  the  sound 
speed  profiles  between  200m  and  2450m.  The  routine  requires  evenly 
spaced  data  and  an  even  number  of  input  points.  The  data  are  there¬ 
fore  interpolated  using  a  cubic  spline  every  50m  creating  46  data 
points . 

Representative  profiles  were  fit  using  the  polynomial  and  the  RMS 
was  computed  for  each  degree  of  fit.  The  objective  was  to  have  the 
mean  RMS  no  larger  than  1  m/sec.  The  minimum  degree  of  fit  meeting 
that  requirement  is  seven  producing  eight  coefficients. 

The  same  technique  is  used  for  salinity  from  0-380m  with  a  three 
degree  fit. 

11  le  deep  portion  of  the  sound  speed  profile,  2500m  to  the  bottom, 
is  very  smooth  and  a  non-orthogonal  least  squares  parabolic  fit  with 
three  coefficients  suffices.  This  fi?  is  actually  made  starting  at 
2000m  to  provide  overlapping. 
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The  near  surface  (0-400m)  temperature  is  the  most  complicated 
with  typical  profiles  simulating  a  step  function.  For  this  fit  an 
analytical  expression  for  the  squared  amplitude  response  of  a  Butter- 
worth  filter  was  developed.  This  technique  works  quite  well  with  most 
profiles  being  fit  with  an  RMS  of  less  than  0.3°C  and  no  rippling  below 
the  mixed  layer. 

SPATIAL  INTERPOLATOR 

A  three-stage  interpolator  is  used  to  produce  matrices  of  the 
coefficients.  The  first  stage  assigns  each  value  to  a  grid  point 
weighted  as  a  function  of  the  inverse  square  of  the  distance  from  the 
grid  point.  If  any  data  fall  within  0.1  grid  interval  of  the  grid 
point  those  data  are  arithmetically  averaged  and  others  excluded. 

The  second  stage  builds  a  coarse  smooth  grid  at  3  times  the  re¬ 
quested  grid  spacing.  This  uses  a  minimum  curvature  cubic  spline 
technique  which  fits  a  surface  to  the  data.  The  technique  was  de¬ 
veloped  by  I.  C.  Briggs  (1974)  and  programmed  by  C.  J.  Swain  (1976). 

The  coarse  grid  is  filled  in  using  a  cubic  spline  to  produce  a  smooth 
surface  at  the  final  grid  spacing  of  30  minutes  in  latitude  and 
longitude . 

The  third  stage  merges  this  smooth  grid  with  the  input  data.  The 
resultant  matrix  shows  detail  where  there  are  data  and  is  smooth  where 
data  are  sparse  with  continuity  of  the  first  and  second  partial  deriva¬ 
tives  maintained  throughout. 

COMPRESSION  AND  RETRIEVAL 

After  the  grids  have  been  evaluated  for  bad  points,  they  are  run 
through  a  series  of  compression,  sorting,  and  indexing  routines.  In 
preparation  for  this,  the  grids  are  organized  so  that  each  season  of 
the  near  surface  includes  the  bathymetry  grid,  the  salinity  coefficient 
grids,  and  the  temperature  coefficient  grids  for  that  season.  Thus 
each  season  consists  of  a  complete  set  for  retrieval.  The  deeper  sound 
speed  coefficient  grids  are  also  combined  with  bathymetry.  This  re¬ 
sults  in  six  compressed  and  sorted  files:  four  seasonal  near— surface 
files  and  two  seasonal  deep  files. 
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Retrieval  routines  simply  require  a  position  and  season  and  the 
appropriate  profile  is  returned.  To  get  a  surface  to  bottom  profile 
requires  accessing  both  a  near-surface  file  and  a  deep  file. 

Therefore,  to  simplify  usage,  a  user  file  is  built  for  each  season 
which  consists  of  reconstructed  sound  speed  profiles  from  surface  to 
bottom  and  salinity  and  temperature  profiles  to  400m  at  30'  positions. 
These  exist  in  the  user  file  similar  to  an  oceanographic  station  with 
values  at  standard  depths. 

The  creation  of  a  complete  sound  speed  profile  from  surface  to 
bottom  requires  two  merges — one  at  400m  and  the  other  at  2500m,  con¬ 
necting  the  three  parts  of  the  model.  The  sound  speed  from  the  surface 
to  400m  is  derived  from  the  temperature  and  salinity  models  and  is 
merged  using  a  modified  version  of  the  ICAPS  merge.  This  merge  con¬ 
sists  of  accenting  the  upper  profile  and  shifting  the  deep  to  fit, 
applying  a  correction  which  decreases  with  depth.  In  order  to  account 
for  the  different  time  steps  between  the  surface  model  and  the  deep 
model,  this  merge  has  been  modified  to  apply  the  same  corrections  in 
reverse  to  the  upper  profile  when  the  difference  between  the  two  at 
400m  is  large. 

The  deep  merge  is  accomplished  with  a  10  point  overlap  and  dif¬ 
ferential  weighting.  This  takes  place  between  2000  and  2450m  with  a 
point  every  50m.  At  20U0m  the  upper  profile  is  weighted  10  and  the 
deep  profile  1  while  at  2450m  the  upper  profile  is  weighted  1  and  the 
deep  profile  10. 

REFERENCES 

Briggs,  I.  C.,  "Machine  Contouring  using  Minimum  Curvature"  Geophysics, 
Vol .  39:1,  pp.  39-48,  1974. 

Swain,  C.  J.,  "A  Fortran  IV  Program  for  Interpolating  Irregularly 
Spaced  Data  using  the  Difference  Equations  for  Minimum  Curvature" 
Computers  and  Geosciences,  Vol.  I,  pp.  231-240,  1976. 
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Attachment  2 

INITIAL  COMPARISONS  BETWEEN  GDEM  AND  TYPICAL  OBSERVED  PROFILES* 

This  attachment  presents  the  preliminary  results  of  qualitative 
tests  to  examine  how  well  the  Generalized  Digital  Environmental  Model 
(GDEM)  reproduces  accurate  historical  sound-speed  profiles  for  numer- 
ical  acoustic  modeling.  For  expediency,  preselected  "typical"  sound- 
speed  profiles  from  the  Northeast  Pacific  region  were  compared  with 
the  GDEM  output.  The  typical  profiles  were  obtained  by  a  NOSC  sta¬ 
tistical  procedure  (described  by  Colbom  and  Pugh,  NUC  TN  1006,  May 
1973)  that  selects  an  observed  profile  to  represent  all  data  for  a 
particular  sound-speed  province  over  a  particular  season.  Figure  1 
shows  the  six  Northeast  Pacific  provinces  (delineated  by  dotted  lines) 
and  the  locations  of  the  actual  profiles  chosen  to  represent  them 
(black  dots  indicated  by  arrows)  over  the  summer  months  of  July  through 
September. 

Each  typical  profile  location  was  given  to  NAVOCEANO  as  input  to 
GDEM.  Figures  2  through  7  plot  the  model's  output  of  summer  profiles 
for  the  six  locations  with  the  real  typical  profiles  from  the  same 
locations.  Quantitative  sound-speed  differences  at  standard  depths 
are  tabulated  on  each  figure.  We  gauged  acceptable  variability  by 
comparing  the  tabulated  differences  with  the  standard  deviation  for 
the  province  data  set  at  each  standard  depth  (standard  deviation  values 
are  not  shown  here).  The  comparison  indicated  good  fittings  for  prov¬ 
inces  2  and  2T.  In  the  other  provinces,  the  typical-GDEM  difference 
exceeds  the  standard  deviation  for  all  province  data  over  the  entire 
season  at  a  number  of  depths. 

A 

This  attachment  is  extracted  from  "Preliminary  Evaluation  of 
Candidates  for  the  Standard  Ocean  Retrieval  System  and  Next  Steps 
Toward  Implementation,"  memorandum  from  the  present  authors  to  Lcdr. 
Kirk  Evans  and  John  H.  Locklin,  JO  April  1979. 

1* 

These  tests  are  limited  to  qualitative  evaluations  because  the 
real  data  for  sound-speed  variability  at  each  location  are  unavailable. 
This  problem  is  being  addressed  in  additional  testing  now  under  way. 
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The  profile  mismatches  are  greatest  at  depths  of  250-500  m  in 
provinces  1,  3,  4,  and  5.  In  provinces  1,  3,  and  5  the  model  profiles, 
though  divergent  from  the  typical  profiles,  still  retain  the  basic 
vertical  structure  of  the  observed  data.  In  province  4  the  problem 
is  more  serious  because  the  model  profile  shifts  the  main  channel 
axis  depth  from  300  m  to  500  m.  The  axis  depth  statistics  for  the 
140  observed  summer  profiles  in  province  4  indicate  a  mean  depth  of 
253  m  and  a  standard  deviation  of  109  m.  Only  8  of  the  140  profiles 
indicate  a  channel  axis  depth  of  500  m,  and  all  occur  near  the  southern 
boundary  of  the  province,  45°-46#N. 

The  deviation  of  the  GDEM  profile  from  the  observed  profile  is 
vividly  displayed  in  Fig.  8.  It  overplots  the  GDEM  profile  on  the 
composite  plot  of  al.  140  observed  summer  profiles  in  province  4. 

1310  region  where  the  surface  (0-400  m)  and  deep  (200-2450  m)  GDEM 
sound-speed  models  overlap  is  shown  for  reference.  The  GDEM  profile 
gradient  between  300  and  500  m  deviates  from  the  majority  of  the  ob¬ 
servations  and  results  in  the  deeper-than-expected  channel  axis  depth. 

Tiie  problem  range,  250-500  m,  contains  the  most  depth  points 
where  the  upper  Butterworth  filter  model  of  the  data  merges  with  the 
middeprh  orthogonal  polynomial  model.  The  nature  and  extent  of  this 
potential  problem  with  CDEM  will  require  additional  testing;  these 
comparisons  merely  suggest  that  GDEM  should  be  evaluated  further  to 
determine  the  quality  of  its  output  for  LRAPP  acoustic  modeling. 


Northeast  Pacific,  area  4,  satraer 
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Attachment  3 

INTERIM  RESULTS  OF  GDEM  ■■  -T-N 
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BAY  ST.  LOUIS,  MISSISSIPPI  39522 

Code  3300: smc 
27  June  1979 


MEMORANDUM 

From:  NAVOCF.ANO  (Code  3300/T.Dav is) 

To:  NOSC  (Joe  Co  I  horn) 

Sub}:  Results  t  rom  new  temperature  model  (0-400M) 

Enel:  (l)  Temperature  comparison  for  PAC  Area  4-Summer 

( 2 )  Sound  Speed  Comparison 

!.  The  summary  you  sent  me  of  our  17  May  meeting  Is  fine.  Based  on 
recent  testing  of  a  new  temperature  model  containing  a  replacement  of 
the  quadratic  tail  by  an  exponential  form,  we  have  "sade  a  major  change 
in  >>ur  1/  May  plan.  instead  of  proceeding  directly  to  the  Atlantic 
model  we  are  rebuilding  the  entire  Pacific  sound  speed  model  with  this 
new  temperature  model.  This  should  be  completed  in  6-8  weeks  and  be 
ready  tor  your  evaluation  work. 

2.  Enclosure  (l)  compares  the  old  and  new  model  using  the  temperature 
data  you  sent  me  tor  SEP AC  area-4  (summer).  Enclosure  (2)  Is  the 
equivalent  sound  speed  using  vour  salinities,  including  a  plot  of 
our  mi  idle  model  sound  speed  before  the  merge  process.  You  can  see 
t ha :  t he  problem  with  the  merge  which  caused  the  middle  model  to  be 
deformed  to  au t c h  file  tall  of  the  surface  model  was  actually  caused 
by  t  tie  poor  f  it  of  the  quadratic  tali.  Test  results  to  date  indicate 
that  t  fie  exponential  tail  on  the  new  temperature  model  will  greatly 
Improve  t  fie  merge  problem. 

>.  Tfier e  is  one  additional  change  to  our  I’  May  plan.  You  and  John 
i.ockliti  convinced  me  ot  t fie  advantages  of  a  separate  temperature  and 
salinity  model  tor  t fie  middle  model.  We  don't  have  time  now  to 
ho  i  !  J  fh.-se  r  or  the  pacific  but  ve  plan  to  take  this  approach  in 
the  At  1  ant  lc . 

•  Thanks  for  sending  me  the  Pacific  listings  and  I'll  let  you  know 
as  soon  as  we  have  the  new  ?a*'iflc  model  finished. 

/'£ — 

TUN  DAVIS  ' 


Lopv  to; 

LEDS  E.E.  Evans  (NOStDA/Codc  600) 

J .  toe  kiln.  OUS l 


□  Input  tamp  points 
o  Old  modsi  (quadratic  taM) 

—  New  modal  (exponential  taM) 


Depth  (m) 
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Sound  speed  (m/sec) 


Enel  (2) 
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Appendix  D 
GFDL* 


GFDL  was  developed  by  the  Geophysical  Fluid  Dynamics  Laboratory 
at  Princeton  University.  The  model's  objective  analysis  of  oceano¬ 
graphic  data  derives  from  the  iterative  difference-correction  method 

••  f 

developed  by  Bergthorsson  and  Doos  and  modified  by  Cressman.  GDFL 

was  designed  to  study  the  ocean's  role  in  the  global  heat  balance. 

Its  output  products  are  analyzed  horizontal  fields  of  temperature, 

salinity,  sigma-T,  and  oxygen.  The  products  are  given  at  all  standard 

depths  between  the  surface  and  5000  m  on  a  l°-square  grid  for  each 

month  over  all  oceans.  It  is  possible  to  compute  sound-speed  values 

and  their  vertical  profiles  at  all  grid  points  down  to  5000  m  from 

the  output  products.  Because  of  its  products  and  worldwide  coverage, 

GFDL  qualifies  as  a  candidate  for  Standard  Ocean. 

ANALYSIS+t 

The  value  for  each  1°  square  is  defined  as  a  mean  representative 
of  the  square's  center.  The  360  x  180  grid  points  are  located  at  the 
intersection  of  1/2°  lines  of  latitude  and  longitude.  In  the  analysis, 
the  average  distance  between  data  points  (the  observed  means)  is 


This  appendix  is  based  on  information  provided  in  a  letter  to 
Daubin  from  Colborn,  29  March  1979,  and  on  the  following  publications: 
P.  Bergthorsson  and  B.  Doos,  "Numerical  Weather  Map  Analysis,"  Tellus, 
Vol.  7,  No.  3,  1955,  pp.  329-340;  G.  P.  Cressman,  "An  Operational  Ob¬ 
jective  Analysis  Scheme,"  Mon,  Wea.  Rev.,  Vol.  87,  No.  10,  1959,  pp. 
367-374;  and  S.  Levitus  and  A.  H.  Oort,  "Global  Analysis  of  Oceano¬ 
graphic  Data,"  Bull.  Amer.  Met.  Soe.,  Vol.  158,  No.  12,  1977,  pp.  1270- 
1284. 

f 

See  the  citations  in  the  preceding  note. 

'ft 'ft 

N0DC  standard  depths  of  0,  10,  20,  30,  50,  75,  100,  125,  150, 
200,  250,  300,  400,  ....  1500,  1750,  2000,  2500,  3000,  4000,  and 
5000  ra.  In  addition,  the  GFDL  analysis  was  performed  at  3500  and 
4500  m. 

++ 

This  description  is  paraphrased  with  permission  from  Levitus 
and  Oort,  who  describe  the  GFDL  analysis  and  its  results  in  detail. 
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computed  by  taking  the  square  root  of  the  total  ocean  area  divided 
by  the  number  of  data  points.  An  influence  radius  is  defined  as  a 
multiple  of  the  average  distance  between  observation  points.  At  the 
grid  points  coinciding  with  an  observed  mean,  the  difference  between 
the  mean  and  the  first  guess  field  is  computed.  At  each  such  grid 
point,  a  correction  to  the  first-guess  value  is  then  computed  as  a 
distance-weighted  mean  of  all  first-guess  mean  differences  occurring 
within  the  grid  point's  influence  radius.  Mathematically,  the  correc¬ 
tion  factor  is  given  by  the  expression 


s=l 


W  Q 
s  s 
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the  correction  factor  at  point  (i,j), 

east-west  and  north-south  coordinates,  respectively,  of 
a  grid  point, 

the  number  of  observations  within  the  influence  radius 
of  point  i,j, 

the  difference  between  the  first  guess  and  the  observed 

mean  at  the  sth  point  in  the  influence  area, 

2  -2 

exp(-Er  R  ),  for  r  £  R, 

0,  for  r  >  R, 

distance  of  the  observation  from  the  grid  point, 
influence  radius, 

4. 


At  each  grid  point,  an  analyzed  value  G  is  computed  as  the  sum 

^  »  J 

of  the  first  guess  F  and  the  correction  C,  ,  at  the  point.  As 

i » J  i  » J 

another  gross  error  check,  if  the  magnitude  of  exceeds  a  prescribed 
limit,  held  constant  throughout  the  analysis,  the  correction  factor 
is  not  used.  If  there  arc  no  data  points  within  the  influence  radius, 
the  correction  is  zero  and  the  analyzed  value  is  simply  the  first-guess 
value . 


Following  the  foregoing  procedure,  an  analyzed  field  is  produced 
at  each  grid  point.  That  field  becomes  the  new  first-guess  field,  and 
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the  procedure  is  repeated  until  the  analysis  shows  little  change. 

The  number  of  iterations  required  depends  on  the  disparity  between 
the  first-guess  field  and  the  observed  input  data.  The  smaller  the 
difference,  the  fewer  iterations  are  required.  After  each  iteration, 
the  resulting  field  is  smoothed  with  a  Laplacian  smoother. 

The  analysis  scheme  permits  the  influence  radius,  the  error  limit 
for  the  correction  factor,  the  number  of  smoothings,  and  the  intensity 
of  smoothing  to  be  varied  with  each  iteration.  The  strategy  is  to 
begin  the  analysis  with  a  large  influence  radius  and  decrease  it  with 
each  iteration.  Levitus  and  Oort  were  thus  able  to  analyze  progres¬ 
sively  smaller-scale  phenomena  with  each  iteration. 

EVALUATION 

GFDL's  basic  data  sources  for  the  original  analysis  are  outdated. 
The  analysis  of  temperature  was  made  with  1.2  million  soundings  com¬ 
posed  of  Nansen  cast  data,  mechanical  bathythermograph  (MBT)  data,  and 
expendable  bathythermograph  (XBT)  data,  all  from  a  pre-1973  NODC  file. 

A  new  analysis  is  planned  to  include  all  Nansen  cast /salinity  temper¬ 
ature  depth  recorder  (STD),  XBT,  and  MBT  data  available  from  NODC  as 
of  June  1976. 

For  Standard  Ocean  application,  GFDL's  most  serious  flaw  is  the 
lack  of  attempt  to  preserve  the  vertical  integrity  in  the  parameter 
fields.  The  iterative  difference-correction  method  is  applied  exclu¬ 
sively  at  horizontal  levels  to  produce  highly  smoothed  horizontal 
parameter  surfaces.  It  ignores  the  unequal  spatial  and  temporal  dis¬ 
tributions  of  data  at  successive  levels.  The  use  of  recombined  param¬ 
eters  at  any  grid  point  to  recreate  vertical  structure  distorts  the 
gradients  and  obscures  acoustically  significant  features,  such  as  sonic 
layer  depth  and  secondary  channels,  in  the  resulting  sound-speed 
profiles. 

an  additional  problem,  recognized  by  Levitus  and  Oort,  is  that 
seasonal  aliasing  is  produced  in  squares  with  incomplete  seasonal  data. 
Given  the  temporal  data  distribution  available,  aliasing  would  greatly 
distort  the  monthly  resolution  in  the  shallow  layers.  Even  with  three- 
month  seasonal  resolution,  the  structure  in  the  final  model  would  be 
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questionable  in  some  regions.  Finally,  the  analvsis  extends  only 
to  a  depth  of  5000  m.  (This  is  not  considered  to  be  a  serious  draw¬ 
back,  however.) 

AVAILABILITY 

The  model  is  presumably  completed  and  available,  although  the 
status  of  the  new  analysis  with  the  revised  data  base  is  unknown. 

The  computer  code  for  the  objective  analysis  was  specifically  developed 
for  the  GFDL  computer,  and  the  output  data  format  may  also  be  computer- 
specific.  The  data  could  be  reformatted,  but  future  updates  of  the 
model  could  not  be  produced  independently  if  it  is  impractical  to  con¬ 
vert  the  code.  Time  and  cost  to  acquire  the  model  cannot  be  estimated 
without  further  investigation.  Computer  compatibility  and  the  labor¬ 
atory's  cooperativeness  are  unknown  factors. 

SUMMARY 

The  GFDL  model  is  not  a  suitable  candidate  for  Standard  Ocean. 

The  objective  analysis  makes  no  attempt  to  preserve  the  integrity  of 
the  vertical  structure  of  the  sound-sneed  field.  Computer-related 
difficulties  prohibit  updating  the  model  in  the  future.  GFDL’s  pri¬ 
mary  advantages  are  that  the  initial  analysis  is  complete,  the  cover¬ 
age  is  good,  and  the  model  is  presumably  available  now. 
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Appendix  E 
HYDAT* 


HYDAT  (Hydro-Climatological  Data  Base)  is  a  software  system  de¬ 
signed  to  retrieve  the  most-representative  observed  surface-to-bottom 
temperature  and  salinity  profiles  for  a  particular  location  and  time, 
based  on  a  data  set  within  a  specified  radius  of  influence  of  the 

location.  Originally  developed  by  Ocean  Data  Systems,  Inc.,  the  pro- 

+ 

gram  has  undergone  much  modification  and  debugging.  According  to 
Evelyn  A.  Hess,  its  operation  is  not  yet  routine  nor  is  HYDAT  ready 
for  release. 

HYDAT  is  intended  to  provide  worldwide  coverage  over  all  months 
for  which  data  are  available.  Hess  believes  it  to  be  the  best  program 
for  retrieving  data  on  the  southern  hemisphere.  For  worldwide  cover¬ 
age,  Hess  considers  the  present  data  set,  an  older  one  of  about  750,000 

observations,  inadequate.  The  installation  of  the  MOODS  data  set  at 

** 

FNWC  in  late  1979  should  improve  the  situation.  MOODS  will  contain 
«  2.8  x  10 k  observations  worldwide. 

HYDAT' s  analysis  relies  on  a  variable,  user-specified  radius  of 
influence.  The  area  analyzed  is  in  1°  x  1°  increments;  smaller  grid 
increments  are  not  used.  Time  spacing  is  one  or  two  months,  though 
smaller  spacings  are  possible  given  adequate  data  density.  Retrieval 
time  is  said  to  be  about  one  minute  per  retrieval,  depending  on  the 
size  of  the  data  set  analyzed. 

The  program  selects  the  mean  surface-to-bottom  temperature  pro¬ 
file  according  to  three  criteria: 

* 

This  appendix  draws  on  discussions  between  Evelyn  A.  Hess 
(FNWC),  Colborn,  and  Daub  in  held  15  March  19  79  at  FNWC,  Monterey, 
California. 

^Ocean  Data  Systems,  Inc.,  Iliulroalima  to  logical  Data  Retreival 
Program:  Functional  Description,  1  November  1976. 

** 

According  to  Lcdr.  Will  Rogers,  FNWC,  the  data  files  for  MOODS 
have  been  assembled,  and  MOODS  is  expected  to  be  on-line  shortly. 
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1.  The  profile  closest  to  the  median  or  mean  sea  surface 
temperature. 

2.  The  profile  closest  to  the  median  or  mean  temperature  at 
243  m  (797  ft). 

3.  The  profile  closest  to  the  mean  heat  content  (average  tem¬ 
perature)  between  0  and  243  m. 

HYDAT's  analysis  ignores  sonic  layer  depth,  which  may  therefore 
not  be  preserved.  Furthermore,  there  is  no  attempt  to  preserve  the 
tiorizontal  integrity  in  the  outputs  from  adjacent  locations.  If  there 
are  fewer  than  six  profiles  within  the  radius  of  influence,  HYDAT 
presents  all  the  data  and  lets  the  user  decide  among  them.  HYDAT  can 
distinguish  and  identify  two  water  masses  within  the  same  radius  of 
influence,  unless  their  temperature  profiles  cross  or  converge  at 
depth.  But  more  than  two  water  masses  are  too  much  for  HYDAT's  analy¬ 
sis.  The  subroutine  Single-SAL  selects  the  most  typical  salinity  pro¬ 
file.  in  the  older  climatology  salinity  data  are  less  abundant  than 
temperature  data. 

SUMMARY 

HYDAT  is  an  unsuitable  Standard  Ocean  candidate  for  the  following 
reasons : 

1.  The  analysis  may  fail  to  preserve  sonic  layer  depth  and 
horizontal  continuity  with  adjacent  areas. 

2.  The  code  is  not  yet  ready  for  release;  an  experienced 
operator  is  required  to  run  the  model. 

3.  The  model  relies  on  an  older,  less  reliable  climatology. 

4.  The  model  may  not  retrieve  representative  profiles  in  areas 
or  during  periods  of  sparse  data. 


-83- 


Append  lx  F 
ICAPS* 


The  Integrated  Command  Antisubmarine  Warfare  Prediction  System 
(ICAPS)  is  an  operational  oceanographic  data  analysis  system  developed 
by  the  Naval  Oceanographic  Office  (NAVOCEANO).  It  is  being  used  by 
units  in  the  fleet  for  sonar  range  prediction.  The  system  can  be  run 
on  a  Univac  1108,  Nova  800,  or  IBM  360  computer. 

ICAPS  merges  observed  XBT  profiles  with  the  appropriate  deep  tem¬ 
perature  profiles  through  a  temperature/salinity  analysis  using  his¬ 
toric  salinities  (water-mass  identi i ication) .  Salinities  above  400  m 
can  be  adjusted  to  maintain  a  nonnegative  local  density  gradient  in 
areas  affected  by  temperature  inversions.  The  operator  can  override 
the  system  if  he  disagrees  with  a  decision.  The  system  contains  a 
library  of  typical  observed  XBT  casts  and  known  seasonal  water-mass 
characteristics  for  the  northern  hemisphere  and  the  Indian  Ocean  to 
20  °S. 

Using  the  merged  temperature  and  historic  salinity  profiles, 

ICAPS  generates  the  local  sound-speed  profile  from  the  surface  to  the 
bottom.  Then  it  computes  the  acoustic  eigen-rays  and  propagation 
losses  in  the  Immediate  vicinity  using  tne  FACT  algorithm.  In  the 
latter  computations  it  is  assumed  that  the  water  mass  is  uniform  near 
the  XBT  cast.  ICAPS  cannot  yet  predict  the  range-dependent  sound- 
speed  structure  from  a  given  point,  but  W.  R.  Floyd  expects  a  predic¬ 
tive  capability  to  be  developed  sometime  in  the  future. 

* 

This  appendix  draws  on  discussions  between  William  R.  Floyd 
(NAVOCEANO),  Paul  M.ersdorf  (NAVOCEANO) ,  Colborn,  and  Daub  in  held  8 
March  1979  at  NAVOCEANO,  Bay  St.  Louis,  Mississippi,  and  on  the  follow¬ 
ing  publications:  Alvan  Fisher,  Jr.,  The  ICAPS  Water  Mass  Histom  File, 
NAVOCEANO,  NOORP-19,  May  1978,  and  idem,  Oceanographic  Analysis  Sianual 
for  On-Scene  Preiliction  Systems,  NAVOCEANO,  N00RP-20 ,  May  1978. 
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METHODOLOGY 

The  following  paragraphs  describe  how  ICAPS  selects  the  water- 

* 

mass  characteristics  and  merges  the  temperature  profiles. 

Two  assumptions  were  made  in  developing  the  ICAPS  historic  water- 
mass  file:  (1)  near-surface  water  masses  can  be  uniquely  identified 
by  thermohaline  characteristics,  and  (2)  the  thermal  characteristics 
of  neighboring  water  masses  are  different  enough  to  permit  reliable 
identification  from  an  expendable  bathythermograph  (XBT)  trace  alone. 
After  identification  of  the  applicable  deep  history,  the  temperature 
values  of  the  input  trace  are  merged  with  deep  temperatures  using  an 
equation  of  the  form 


Tt  “  THt  +  K^K^AT)  ,  (F.l) 

where  and  TH^  are,  respectively,  estimated  and  historical  tempera¬ 
tures  at  depth  i ;  K  is  a  weighting  factor;  and  AT  is  the  difference 
between  the  temperature  at  the  bottou  of  the  XBT  trace  and  the  inter¬ 
polated  historical  temperature  at  the  same  depth.  The  weighting 
factor, developed  from  empirical  solution  for  a  sat  of  historical 
data,  is  determined  as  a  function  of  the  depth  increment  between 

points  (D.  -  D .  , ) : 

t  t  - 1 

(D  -D  . ) / 100 

-  0.815  1  .  (F. 2) 

At  the  first  synthesized  temperature  value  (i  »  1),  equals  unity. 

Because  classical  oceanographic  literature  provides  few  guidelines 
for  water-mass  identification,  it  was  decided  that  the  most  objective 
way  of  determining  water-mass  characteristics  within  a  given  area  was 
to  review  original  oceanographic  data.  Two  NAVOCF.ANO  data  files  were 

A 

This  subsection  is  paraphrased  from  Fisher,  The  ICAPS  U\iU;r  Ikuui 
Hietov'.i  Pile,  pp.  1-iO,  with  permission. 

Later  evaluation  of  the  merge  showed  that  a  constant  of  0,790 
created  a  more  realistic  merge  in  the  Mediterranean  Sea.  The  value 
of  0.835  was  retained  for  all  other  areas. 
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available.  (1)  An  oceanographic  station  data  file  of  approximately 
491,000  observations  compiled  by  the  National  Oceanographic  Data  Center 
(NODC)  provided  temp •  ature  and  salinity  data  at  32  standard  depths 
between  the  surface  and  7000  m.  (2)  An  XBT  file  of  approximately 
218,000  observations  from  three  sources  (NAVOCEANO,  NODC,  and  FNWC) 
provided  temperature  data  at  each  flexure  point  over  the  depth  range 
of  the  instrument  (as  deep  as  760  m) .  The  following  procedure  was  used 
to  determine  water-ma  *.  characteristics  in  the  near-surface  layer 
(0-400  m): 

1.  The  classical  literature  was  searched  for  applicable  de¬ 
scriptions.  For  example,  the  northern  edge  of  the  Gulf  Stream 
is  frequently  delineated  by  the  15°C  isotherm  at  200  ra. 

2.  The  ocean  station  data  file  was  used  to  provide  annual  com¬ 
posite  statistics  (mean,  standard  deviation,  number  of  obser¬ 
vations)  at  each  standard  depth  using  all  available  data 
within  the  area  of  interest.  Plots  of  the  distribution  of 
temperature  versus  salinity  at  200  and  400  m  helped  determine 
the  number  of  water  masses  present  and  the  thennohaline  vari¬ 
ability  within  each.  Figure  F.l  shows  a  plot  of  temperature 
versus  salinity  at  200  m  in  an  area  where  the  cold  Labrador 
current  meets  the  warmer  North  Atlantic  drift.  The  presence 
of  water  masses  with  specific  thermohal ine  characteristics 

is  clearly  recognizable,  and  tentative  water-mass  classifica¬ 
tion  has  been  made.  The  200  m  level  was  found  to  be  a  good 
depth  for  classification  since  it  is  too  deep  for  diurnal  and 
seasonal  influences  ytt  within  the  depth  range  of  XBT  probes. 
Tiie  XBT  file  provided  statistical  data  and  histograms  for 
temperature  and  temperature  gradients  at  preselected  depths 
Co  supplement  the  ocean  station  data. 

).  Flexure  points  in  the  temperature  versus  salinity  (T-B)  plot 
shown  In  Fig.  F.l  clearly  defined  water-mass  criteria  in 
areas  where  different  water  masses  existed  in  close  proximity. 
Considerable  temperature  variability  also  occurred  in  areas 
containing  a  single  water  mass,  probably  a  result  of  dynamic 
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Fisher,  The  !  'At':'  Hitter  Mm«  HioU>r[,  Fite,  p.  3  (used 
by  permission). 


Fig.  F.l--0istr1bution  of  temperature  versus  salinity 


-87- 


events  such  as  upvelling.  Where  variability  of  this  nature 
was  observed,  two  classifications  ("warn"  and  "cold")  were 
made  to  provide  a  better  merge  between  XBT  trace  and  history. 

4.  Temperature  ranges  (filters)  at  200  m  were  developed  to  dis¬ 
tinguish  adjacent  water  masses  based  on  information  provided 
in  the  previous  steps.  If  adjacent  water  masses  had  similar 
temperature  ranges  at  200  m,  they  were  differentiated  by 
examining  the  temperature  gradient  between  200  and  300  m. 

For  example,  both  the  Gulf  Stream  and  the  Sargasso  Sea  are 
characterized  by  a  temperature  range  of  15*  to  25*C  at  200  m. 
Analysis  of  a  near-isothermal  layer  of  18*C  water  extending 
from  the  bottom  of  the  seasonal  thermocline  to  over  300  m  deep 
in  the  Sargasso  water  far  from  the  Gulf  Stream  showed  that 

95  percent  of  the  observations  had  a  temperature  gradient 
between  0.0*C/100  m  and  -1.6*C/100  m.  Thus,  in  the  region 
of  the  Gulf  Stream,  the  gradient  -1.6*C/100  u  at  the  200- 
300  m  level  is  used  to  differentiate  Gulf  Stream  water  from 
Sargasso  water. 

5.  Mean  seasonal  temperature  and  salinity  values  were  then  de¬ 
termined  for  each  depth  and  water  mass  (Table  F.l).  Where 
the  data  were  not  deep  enough,  temperature  and  salinity  were 
extrapolated  to  the  bottom  by  comparison  with  neighboring 
profiles.  Inconsistencies  In  the  data — sueh  as  a  temperature 
inversion  at  depths  below  200  a — were  examined  to  determine 
if  thev  were  a  result  of  statistical  processing,  data  dis¬ 
tribution,  or  bad  data. 

6.  A  quality  control  cheek  was  made  by  plotting  the  seasonal 
data  on  a  single  plot  of  temperature  versus  salinity  (Fig- 
F.2).  That  procedure  immediately  reveals  inconsistencies  In 
the  data:  temperature  errors  are  indicated  by  vertical  spikes, 
salinity  errors  by  horizontal  spikes,  and  depth  errors  by 
skewed  spikes.  Where  data  were  obviously  incorrect,  the  plot 
was  smoothed  to  conform  with  surrounding  data.  A  second 
quality  control  check  was  made  by  visually  inspecting  the 
seasonal  traces  of  temperature  and  salinity  versus  deotb. 

Again,  discrepancies  were  smoothed  after  eoc^parison  with 
neighboring  traces. 
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Table  F.  1 

TEMPERATURE  AND  SALINITY  AT  STANDARD  DEPTHS  IN  SLOPE  WATER 


Standard 

Depth 

Temperature  (*C) 

Salinity  ( 

*/..) 

Mean 

Standard 
Devlat ion 

Numher 

of 

Observat ions 

Mean 

Standard 

Deviation 

Number 

of 

Observations 

0 

23.77 

2.28 

676 

34.34 

1.10 

684 

10 

23.10 

2.65 

682 

34.55 

.97 

680 

20 

21.72 

3.70 

682 

34.84 

.83 

680 

30 

19.42 

4.36 

683 

34.96 

.83 

679 

50 

15.87 

4.17 

683 

35.14 

.76 

678 

73 

14.49 

2.86 

683 

35.41 

.55 

678 

100 

13.78 

2.04 

683 

35.54 

.38 

678 

125 

13.16 

1.57 

684 

35.54 

.28 

678 

i50 

12.54 

1.34 

684 

35.51 

.22 

678 

200 

11.21 

1.24 

684 

35.38 

.17 

676 

250 

9.86 

1.23 

684 

35.24 

.15 

674 

300 

8.68 

1.25 

682 

35.14 

.13 

674 

400  | 

6.87 

i.10 

582 

35.03 

.10 

575 

500 

5.67 

.79 

551 

34.99 

.06 

546 

600 

5.03 

.52 

529 

34.98 

.05 

525 

700 

4.67 

.  32 

518 

34.98 

.04 

514 

800  | 

4.43 

.24 

473 

34.97 

.03 

471 

900 

4.27 

.20 

4  38 

34.97 

.03 

436 

1000 

4.13 

.17 

393 

34.96 

.03 

386 

1100 

4.02 

.15 

350 

34.96 

.04 

345 

1200 

3.92 

.13 

330 

34.96 

.04 

324 

1 300 

3.85 

.  12 

322 

34.96 

.04 

316 

1400 

3 .  73 

.12 

319 

34.95 

.04 

313 

l  ‘.00 

3.72 

.12 

315 

34.95 

.04 

311 

P50 

3.  56 

.09 

2  70 

34.  >5 

.04 

264 

20(H) 

3.41 

.09 

2  39 

34.95 

.04 

2  33 

2  300 

3.00 

.11 

160 

34.94 

.03 

154 

3000 

2.59 

.16 

89 

34.92 

.03 

84 

4000 

1 

2.26 

.07 

41 

34.90 

.02 

38 

SOURCE:  Adapted  from  Fisher,  The  i CAST  Water  foists  Hinioty  File,  p.  5 
(used  by  pcrelsaion). 


NOTE:  Slope  water  character  1st  less  are  as  follows:  location,  35*- 

42*N,  60*-?6"U;  season,  summer;  temperature  »age ,  9*-15*C;  aallnl.y 
range,  30-40  */... 
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The  water-mass  file  has  been  segmented  to  permit  installation  in 
computers  of  various  storage  capacity.  For  example,  the  North  Atlantic 
is  divided  into  areas  A  though  E,  the  North  Pacific  A  through  G,  and 
the  Indian  Ocean  A  through  D.  Each  area  is  further  divided  into  regions 
of  similar  oceanographic  properties,  the  lowest  denominator  being  a  1° 
rectangle.  A  region  may  have  as  many  as  five  water  masses  but  is 
normally  limited  to  two  or  three.  Historical  data  are  provided  by 
season.  In  most  regions  winter  is  January  through  March.  In  the  Indian 
Ocean,  however,  the  winter  monsoon  is  October  through  March;  summer 
monsoon,  April  through  September.  Given  the  geographic  position,  data 
from  an  XBT  trace,  and  season,  the  program  will  automatically  select 
the  proper  water-mass  history  for  the  merge. 

EVAI  UAH  ON 

For  several  reasons,  ICAPS  is  an  attractive  candidate  for  Standard 
Ocean: 

1.  The  system  is  highly  automated  and  can  be  operated  success¬ 
fully  by  a  nonoceanographer. 

2.  It  is  operational  and  available. 

3.  It  can  differentiate  and  Identify  water  masses. 

W.  R.  Floyd  considers  ICAFS  a  rel.able  system  in  the  North  Atlantic, 
less  so  in  the  North  Pacific. 

Other  features  of  ICAPS  make  it  unsuitable  for  Standard  Ocean: 

1.  ICAPS  cannot  make  range-dependent  predictions  of  sound  speed. 

It  assumes  a  uniform  environment  in  the  vicinity  of  the  XBT 
cast  input.  Therefore,  it  will  not  function  well  near  fronts 
or  in  highly  variable  regions. 

2.  The  system  requires  XBT  d3ta  as  Input. 

3.  The  identification  of  water  mass  and  subsequent  extension  of 
the  XBT  profile  to  the  bottom  may  not  be  reliable  in  areas 
where  rhe  deep  climatology  (>400  m)  varies  significantly 
(e.g.,  the  eastern  North  Atlantic  near  the  Mediterranean 


outflow).  Having  made  no  assessment  of  ICAPS'  performance, 
the  Standard  Ocean  Evaluation  Group  is  unsure  of  the  serious¬ 
ness  of  the  foregoing  disadvantage. 


sr„ 
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Appendix  G 
ODSI* ** 


Ocean  Data  Systems,  Inc.  (ODSI)  proposes  to  develop  an  environ¬ 
mental  data  retrieval  system  to  be  used  as  Standard  Ocean.  Captain 
P.  M.  Wolff  is  to  direct  five  analysts  in  the  effort.  Captain  Wolff 
plans  to  base  the  model  on  similar  work,  done  by  ODSI  for  the  Department 
of  Energy.  For  each  oceanic  square,  the  model  will  contain  the  monthly 
most  probable,  minimum  and  maximum  temperature,  salinity,  sound-speed 
profile,  and  sonic  layer  depth.  ("Most  probable"  refers  to  typical 
oceanographic  conditions  in  an  area;  it  is  not  an  average  or  a  mean  of 
local  observations,  which  might  portray  an  environment  that  would  never 
be  observed.)  Profiles  will  extend  from  the  surface  to  the  greatest 
depth  in  the  square.  Squares  will  vary  in  size  depending  on  local 
complexity;  typically,  they  will  be  1*  *  1*.  In  addition  to  the  pro¬ 
files,  ODSI  would  provide  the  following  oceanographic  Information  for 
each  square: 


•  Current  and  frontal  boundaries. 

•  The  presence  of  strong  bottom  currents,  upwelling,  and 
eddies . 

•  Shelf  phenomena  (not  defined)  and  shelf  boundaries. 

•  u  profiles  (minimum,  maximum,  and  most  probable). 

•  basin  sill  depths. 

•  Discrete  profiles  for  two  wat.,i  masses,  if  present  (how 

** 

they  would  be  differentiated  is  unexplained  ). 


* 

This  appendix  draws  on  discussions  between  Capt .  Paul  M.  Wolff 
(ODSI),  Colburn,  Daubin,  and  Hashitaoto  held  16  March  1979  at  Ocean 
Data  Systems,  Inc.,  Monterey,  California. 
f 

P.  M.  Wolft,  Tempeixttui'e  Varialh  l  i  t at  Three  OTRC  Gitee,  Ocean 
Data  Systems,  Inc.,  July  1978;  Ocean  Data  Systems,  Inc.,  and  U.S. 
Department  of  Energy,  OT&C  Thermal  Reeout\?e  Report  for  Western  Gulf 
of  Mexico,  T ID- 7 7949 ,  October  1977;  Idem,  OTSC  Thermal  Re&ouroe  Re¬ 
port  for  Central  Gulf  of  Mexxoo,  T1D-279S1,  n.d. 

** 

Appendixes  E  and  F  describe  algorithms  for  identifying  water 


masses . 
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The  data  would  be  arranged  in  tables,  with  eddies,  currents,  and  the 
like  represented  by  code. 

The  data  to  he  used  in  deriving  the  model  will  come  from  the 
sources  in  ODSl's  extensive  in-house  data  library.  Three  files  will 
he  maintained  for  the  model's  preparation. 


Data  Source 


File 

1. 

Nansen  casts  1 

STDs  > 

CTDs*  ) 

NODC 

File 

2. 

XETs  \ 

NODC 

(50%) 

AXBTs  +  J 

FNWC 

(50Z) 

File 

3. 

MBTs 

NODC 

Flit*  1  will  take  precedence  over  files  2  and  3  when  sufficient  data 
are  available.  File  3  will  be  used  only  in  holidays  where  no  data  from 
files  1  and  2  are  available.  If  data  are  unavailable  for  certain 
areas  from  any  file,  values  will  be  extrapolated  from  adjacent  areas. 
Captain  Wolff  will  avoid  analyzing  adjacent  data  sets  containing  data 
from  different  time  periods.  Rather,  he  will  attempt  to  analyze  what 
might  be  called  temporally  averaged  data  sets.  The  space  scale  will 
be  appropriate  to  the  local  situation.  For  example,  in  areas  near 
frontal  boundaries  such  as  the  Gulf  Stream,  grid  spacing  would  be  much 
finer  (<1®  *  1*)  than  in  the  Sargasso  Sea,  where  little  spatial  and 
temporal  variability  is  expected. 

Boundaries  and  boundary  widths  for  currents  and  fronts  would  be 
used  to  set  boundary  conditions,  which  would  be  moved  to  match  ob¬ 
served  conditions  in  the  model's  synoptie  version.  The  model  would 
preserve  horizontal  and  vertical  gradients. 

EVALUATION 

An  Impress i w  amount  and  variety  of  data  are  to  be  included  in 
the  proposed  OOSt  model — even  exceeding  the  base  requirements  for 

Conduct ivit v/temperature/dcpth  instrument  readings. 

Airborne  expendable  bathythermograph  traces. 
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Standard  Ocean  in  some  respects.  As  ODSI  remains  conceptual  and  un¬ 
tested,  however,  the  Standard  Ocean  Evaluation  Group  has  no  way  of 
Judging  the  quality  of  the  model's  output  products.  Running  an  analy¬ 
sis  of  even  a  small  area  for  preliminary  evaluation  would  entail  sub¬ 
stantial  start-up  costs.  Furthermore,  it  is  doubtful  that  the  model 
can  be  satisfactorily  completed  for  the  estimated  cost  ($180,000)  and 
within  the  estimated  time  (18  months).  Other  candidates  such  ns  GDEM 
and  AUTO-OCEAN,  while  lacking  the  quantity  of  output  products  of  ODSI, 
are  available  for  evaluation  now  and  should  be  less  expensive  to 
acquire.  For  those  reasons,  we  recommend  against  adopting  the  ODSI 
model  for  Standard  Ocean. 
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Appendix  11 
SIMAS* 


The  Sonar  In-Situ  Mode  Assessment  System  (SIMAS)  sound-speed 
data  file  was  developed  hy  the  Naval  Uncervater  Systems  Center  (NUSC), 
New  London.  Its  purpose  is  to  indicate  the  sound-speed  structure  below 
366  o.  extending  an  observed  XBT-derived  sound-speed  profile.  The 
basic  file  contains  annual  sound-speed  profiles  from  366  m  to  the 
bottom  for  70  so-called  homogeneous  areas  in  the  North  Atlantic,  69 
areas  in  the  North  Pacific,  and  50  areas  in  the  Indian  Ocean.  Cover¬ 
age  for  the  Atlantic  and  Pacific  extends  from  63*N  to  10*S;  Indian 
Ocean  coverage  extends  from  the  coast  to  50°S.  The  data  have  been 
digitized  geographically  to  a  presumed  resolution  of  at  least  1*  lati¬ 
tude  and  1°  longitude.  In  addition,  monthly  best-guess  sound-speed 
profiles  are  available  for  the  upper  366  m,  although  they  may  not  be 
in  a  spatially  digitized  format.  The  complete  profiles  are  available 
in  atlas  reference  charts  by  area.  The  program  and  data  are  computer¬ 
ised  for  shipboard  operation.  The  system  resembles  ICAPS  (Appendix 
F)  in  purpose,  except  that  ICAPS  can  identify  the  local  water  mass 
whereas  SIMAS  depends  on  fixed  sound-speed  province  boundaries. 

EVALUATION 

The  SIMAS  analysis  is  based  on  the  original  NODC  Nansen  cast  data 
base  available-  bet'or-.-  1976,  when  th  original  analysis  was  reported. 

The  data  were  edited  by  lteepi*ig  only  the  most  recent  observations 
while  trying  to  retain  spatial  and  temporal  coverage.  Through  a  sub¬ 
jective  analysis,  standard  profiles  and  area  boundaries  were  selected. 
After  analysis  the  data  set  was  reduced  to  only  0.1  percent  of  the 

-  ~~  - 

This  appendix  craws  on  Information  in  a  letter  .o  Daub  in  from 
Colborn,  29  March  1979. 

♦ 

Eugene  M.  Podeszwa,  Sound.  S^-e^d  Profile  for  the  North  pturif to 
'Mean,  Naval  Underwater  Systems  Center,  TD-5271,  2  February  1976; 
idea,  ik>und  Profi  log  for  the  North  Atlantic  Cheitn ,  Naval  Under¬ 

water  Systems  Center,  TD-S447,  20  October  1976;  idea,  Sound  Speed  Pro¬ 
files  for  f'u  Indian  Ansum,  Naval  Underwater  Systems  Center,  TD-S5S5, 

11  Deceober  1976. 


FHECKttllC  P AOS  hUMUitfTT  FILMED 


-98- 


orlglnal  observations.  Vertical  smoothing  of  the  data  was  performed. 
Horizontal  smoothing  was  used  to  produce  a  single  structure  for  the 
entire  North  Pacific  below  2134  m. 

It  is  difficult  to  evaluate  the  technical  quality  of  the  sound- 
speed  profiles  because  of  the  subjectivity  of  the  analysis.  We  can, 
however,  point  out  defects  in  the  homogeneous  sound-speed  area  format 
designed  to  provide  a  single  profile  for  a  given  location  and  month. 
That  structure  is  not  a  true  sound-speed  field;  it  does  not  have 
realistic  horizontal  gradient  characteristics;  and  it  cannot  be  re¬ 
liably  used  to  produce  range-dependent  model  inputs. 

The  availability  of  SIMAS  is  not  known.  However,  the  analysis 

* 

is  complete  and  published  in  atlas  form.  The  time  and  cost  to  acquire 
the  digital  data  file  would  depend  on  its  compatibility  with  the  com¬ 
puter  format  of  the  LRAPP  data  bank.  Acquisition  might  be  accomplished 
in  a  fairly  short  time  and  at  reasonable  cost. 

SUMMARY 

The  SIMAS  data  file  is  not  a  suitable  Standard  Ocean  candidate 
because  It  has  not  been  developed  to  produce  a  sound-speed  field  with 
horizontal  continuity.  Extensive  testing  would  be  needed  to  resolve 
questions  of  Its  technical  quality.  The  model’s  main  advantage  is 
its  presumed  availability  now. 

* . 

See  the  three  Podeszva  reports  cited  above. 
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TO: 

FROM: 


Distribution 


SUBJECT: 


ENCL. 


14  March  1980 


Scott  C.  Daub in,  Jr. 

Coordinator,  Standard 
Ocean  Evaluation  Group 
SEAS  Modeling  Program 

Evaluation  of  Standard  Ocean  Candidates: 
Forwarding  of  Report 

PSR  Report  922,  Evaluation  of  Standard 
Ocean  Candidates ,  March  1980 


1.  The  Long  Range  Acoustic  Propagation  Project  (LRAPP)  of  the  Naval 
Ocean  Research  and  Development  Activity  (NORDA)  has  sponsored  a  search 
f'”-  an  oceanographic  data  retrieval  system,  to  be  installed  as  Standard 
Oc^.an  in  the  LRAPP  data  bank.  (LRAPP  is  now  entitled  the  Surveillance 
Environmental  Acoustic  (SEAS)  Project.)  Standard  Ocean’s  primary  pur¬ 
pose  will  be  to  provide  range-dependent  sound-speed  profiles  for  input 
to  NORDA’s  numerical  acoustic  models.  Standard  Ocean  will  also  support 
the  analysis  of  environmental  data  collected  during  exercises  at  sea. 
Eight  existing  or  proposed  candidate  systems  were  assessed  by  t he 
Standard  Ocean  Evaluation  Group  between  January  and  August  of  1979. 

Each  candidate  was  rated  for  its  ability  to  meet  LRAPP’ s  requirements 
for  accuracy,  ease  of  use,  speed,  availability,  and  cost.  Two  of  the 
eight  candidates,  AUTO-OCEAN  and  GDEM,  were  found  to  meet  or  exceed 
most  Standard  Ocean  criteria. 

2.  The  enclosed  report  by  the  Standard  Ocean  Evaluation  Croup  details 
the  results  of  its  assessment.  The  report  is  authorized  by  the  manager, 
SEAS  modeling  program. 
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Enclosure:  As  noted 

Distribution:  Printed  at  end  of  report 


DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
875  NORTH  RANDOLPH  STREET 
SUITE  1425 

ARLINGTON  VA  22203-1995 

IN  REPLY  REFER  TO; 


5510/1 

Ser  32 1OA/0 11/06 
31  Jan  06 


MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 .  In  accordance  with  reference  (a),  a  declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


^ — fr. 

BRIAN  LINK 
By  direction 


Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
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